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DITRODUCTION  -  LAM31^AR-FL0W  AIRF0II5 


1.  In  high  speed  flight  conditions,  the  drag  of  an  airfoil  is 
almost  exclusively  due  to  skin  friction.   Therefore,  if  further  reduc- 
tion in  drag  is  desired,  it  is  necessary  to  delay  as  much  as  possible 
the  transition  from  laminar  to  turbulent  flow  in  the  boundary  layer  along 
the  surface,  thus  decreasing  the  extent  of  the  turbulent  boundary  layer 
which  gives  considerable  skin  friction.  As  the  factors  that  may  affect 
the  transition,  we  will  consider  the  stream  tiirbulence,  the  surface 
roughness,  the  surface  pressure  distribution,  and  so  on.   In  actual 
flight  conditions,  however,  the  effect  of  turbulence  seems  to  be  unex- 
pectedly small,  so  that,  so  far  as  smooth  surfaces  are  concerned,  there 
remains  only  the  shape  of  the  airfoil  section  in  relation  to  pressure 
distribution  as  the  most  important  factor  affecting  transition.  We  call 
a  laminar-flow  airfoil  that  airfoil  in  which  the  shape  of  the  section 

is  suitably  designed  so  as  to  delay  the  transition  of  the  boundary 
layer  flow. 

2.  It  is  evident  that  the  laminar  separation  of  the  boundary  layer 
may  cause  the  transition,  as  will  be  mentioned  in  the  appendant  part  of 
the  paper,  paragraphs  35-^0.  We  cannot  expect,  therefore,  to  maintain 
laminar  flow  beyond  the  separation  point.   Summarizing  the  results  of 
flight  experiments  on  airfoils  hitherto  made  (refs.  8  to  12),  we  have 
the  conclusion  that  the  observed  transition  coincides  approximately  with 
the  calculated  laminar  separation  point  at  small  Reynolds  n\ambers,  while 
it  moves  upstream  toward  the  minimum  pressure  point  as  the  Reynolds 
number  increases.  However,  no  example  has  ever  yet  been  observed  in 
which  the  transition  moves  ahead  of  the  minimum  pressure  point.   We 
therefore  arrive  at  the  supposition  that  the  laminar-flow  airfoil  may  be 
most  simply  realized  by  designing  the  airfoil  in  which  the  minimum  pres- 
sure occurs  well  downstream. 
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DESIGN  OF  SYMMETRICAL  AIRFOII^  IN  WHICH  THE  MINIMUM 
PRESSURE  OCCURS  DOWNSTREAM 

3.  Following  Professor  Moriya  (ref.  13) >  we  write  the  coordinate 
along  the  chord  in  the  form  x  =  —  (l  +  cos  I),  and  assign  x  =  0, 

I  =  jT  to  the  leading  edge,  and  x  =  1,  1=0  to  the  trailing  edge. 

00 

Expressing  the  ordinate  of  the  mean  camber  line  hy  M  =  2 ^   cos  nl, 

0 

00  ^ 

and  the  half -thickness  measured  normal  to  the  chord  by  T  =  ^   b^  sin  nl  f 

1 
the  pressure  distribution  around  the  airfoil  in  the  two-dimensional 
potential  flow  is  given  by 


=  1 


cos 


I 

a  "j:;  —  sin  I  +  2__  naj^  (l  -  cos  nl)  +  )   nb^  sin  n^ 


+  sm  a  < 


^  (1  -  cos  5)  +  2 ^^n  ^^^  ^^   ^   I ^^   ^°^  ^^ 

11 


1   .  2  t 

r-  sm   ^  +  <| 


+  2 '^^n  sin  nl  +  ^ rib^^  cos  nl 

1  1 


=  1 


cos  a  <|  1  i  fA(,  +  eBg  I  +  sin  a< 


"^ 


1  -  X 


-  \|-T-  -  (^^s  ^  eB,) 


1  +  (fAg  i  eBc)' 


iThe  ordinates  of  the  upper  and  lower  surfaces  are  given  by  M  +  T 
and  M  -  T,  respectively. 
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where  a  "is  the  angle  of  attack,  p  is  the  pressure  acting  on  the  air- 
foil surface,  measured  from  the  static  pressure  of  the  xindisturbed  stream, 
and  q  is  the  dynamic  pressure  of  the  undisturbed  stream.   We  assume  that 

00 

the  trailing  edge  is  sharp,  so  that  2 ^^n   =  *-'•   ^^  limit  the  range  of 

1 
the  variables  I     between  0  and  it,  and  assign  the  upper  and  lower 
parts  of  the  double  sign  for  the  upper  and  lower  surfaces,  respectively. 
Writing  f  for  the  maximum  value  of  M  (the  maximum  camber)  and  e  for 
the  maximum  value  of  2T   (the  maximum  thickness),  we  put 


fAc 


1  -  cos  nl  rr^-  gin  nl     dM 


T 1  -  cos  n^  \  ^ 

-2  t_  na^  -— —       fAg  =  -2  )_  na^ 


sin  I  1      Bin  I      dx 

v^     sin  n|  r-^     cos  n|     dT 

eBg  =  2  2_  nbn  — —  eB^  =  2  ^  nb^  —  =  -  — 

1      sin  5  1      sm  I      "^ 

The  lift  coefficient  is  given  by 


I 

C-]^  =  2rt  J  sin  a  -  2  2 na^^  cos  a 


We  consider  first  only  the  thickness  of  the  airfoil  (the  camber  of 
the  center  line  will  be  considered  in  the  next  section,  paragraphs  8 
to  11) .  Namely,  we  consider  the.  symmetrical  airfoil  section  set  at  zero 
angle  of  attack,  with  a  view  to  obtaining  the  minimum  pressure  well 
downstream. 

k.   We  adopt  as  the  typical  example  of  the  commonly  used  symmetrical 
airfoils  the  NACA  symmetrical  airfoil  (ref .  ik) 


T  =  e 


1.1+845^  -  O.630OX  -  1.7580x2  +  1.1^215x3  -  O.'yoj^x^i 


The  maximum  thickness  is  located  at  x  =  O.3,  the  leading-edge  radius  is 
l.le^,  and  the  trailing-edge  slope  -{dT/dx)^^^.     ^^   l.lTe.  The  pressure 

distribution  for  the  case  e  =  0.1  is  shown  in  figure  1.  The  minimum 
pressure  is  located  at  x  =  0.1,  and  the  laminar  separation  point, 
determined  by  the  approximate  method  due  to  the  author  (refs.  15  and  I6) , 
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at  X  =  0.61.  If  the  transition  point  of  the  boimdary  layer  would  not 
move  upstream  "beyond  the  laminar  separation  point,  we  might  expect  to 
maintain  a  laminar  boundary  layer  for  more  than  half  the  surface  of  the 
airfoil.   The  flight  experiments  hitherto  made,  however,  appear  to  give 
negative  evidence  for  such  a  conjecture. 

5.  Now,  in  order  to  shift  the  minimvrai  pressure  backward,  it  is 
required  to  shift  the  position  of  maximum  thickness  (x  =  m)  backward. 
For  designing  such  airfoils,  we  represent  the  shapes  of  parts  before 
and  after  the  maximum  thickness  by  two  algebraic  expressions.  For  the 

forward  half  (O  ^  x  ^  m) 


/- 


T  =  e 


y2hx  +  h-^x  +   hgX^ 
while  for  the  rear  half  (m  ^  x  =  l) 

T  =  e  <^0.01  +  did  -  x)  +  d2(l  -  x)^  +  d3(l 

where 

2  -  3^2hm  V2hm  -  1 


hi  =  —^ h2  = 


2m'= 


l.i^T  -  2di(l  -  m)  dj^Cl  -  m)  -  O.98 

6.2   =  do  =  ■ 

(1  -  m)2  (1  .  ni)3 

and  we  assign  arbitrary  values  for  three  parameters,  m,   h  (=  leading- 
edge  radiiis  ^  ^  ),   and  di  (=  trailing-edge  slope  -i-  e)  .  Although  the 
method  has  the  drawback  that  the  two  expressions  give  different  values 

of  d'^T/dx^  at  x  =  m,  where  dT/dx  becomes  zero,  we  nevertheless 
adopt  it  because  we  are  in  a  position  to  vary  the  forward  and  rearward 
parts  most  simply  and  independently. 

6.  First,  we  fix  the  forward  half  with  m  =  0.5  and  h  =  0.5,  and 
vary  the  rear  half  by  giving  d-,   the  values  1.7,  2.0,  2.5,  and  3.0, 

respectively.  The  shape  of  the  section  and  the  pressure  distribution 
for  e  =  0.1  are  shown  in  figure  2.  We  find  from  this  result  that,  as 
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di   increases,  the  minimvun  pressure  point  moves  backward  and  the  gradient 

of  pressure  rise  following  the  minimum  pressure  steepens.  We  also  find 
that  the  pressure  distribution  in  the  neighborhood  of  the  minimum  pres- 
sure exhibits  a  wavy  indentation  when  the  value  of  d3_  is  too  small  or 

too  large,  and  that  there  exists  a  certain  value  of  d^^  for  which  the 
pressure  distribution  is  flat  and  smooth.  Such  a  value  of  d^.      ^^   about 
2.5  in  this  case.  We  therefore  fix  the  rear  half  with  d^  =  2.5,  and 

vary  the  forward  half  by  giving  h  the  values  0.35,  O.5O,  O.7O,  and 
1.05,  respectively.  The  shape  of  the  section  and  the  pressure  distri- 
bution for  e  =  0.1  are  shown  in  figure  3-  From  this  comparison,  we 
find  that  the  negative  pressiire  bump  Immediately  behind  the  leading 
edge  decreases  as  h  decreases,  and  that  the  maximiim  permissible  value 
of  h  is  about  0.7. 

The  effect  of  thickness  is  shown  in  figure  k,    in  which  curves  of 
pressiore  distribution  are  given  for  different  values  of  e,  O.O6,  0.10, 
and  O.l^t,  but  for  a  fixed  set  of  parameters,  m  =  O.5,  h  =  0,5,  and 
di  =  2.5.   It  is  seen  that  the  characteristics  of  the  pressure  distri- 
bution do  not  materially  change  with  thickness.  There  is,  however,  a 
slight  change  in  the  pressure  distribution,  the  maximum  permissible 
value  for  h  slightly  increasing  as  the  thickness  increases. 

To  see  the  effect  of  the  position  of  maximum  thickness,  we  give  m 
values  ranging  from  0.35  to  O.6O,  varying  at  the  same  time  values  of  d]_ 
and  h  so  that  the  pressure  distribution  becomes  flat  and  smooth.  The 
result  of  calculation  is  given  in  fig\ire  5,  which  shows  a  considerable 
change  in  the  position  of  minimum  pressure.  The  change  is  not  purely 
due  to  the  effect  of  m,  but  it  is  at  any  rate  to  be  noticed  that  the 
value  of  m  less  than  0.^  is  not  sufficient  for  shifting  backward  the 
minimum  pressure,  while  increasing  the  value  of  m  beyond  0.5  is  of  no 
advantage,  since  the  backward  shift  is  then  almost  saturated,  only  the 
adverse  pressure  gradient  being  increased. 

7.  From  the  results  of  calculation,  we  thus  arrive  at  the  conclusion 
that  m  miist  be  between  O.k   and  O.5  and  h  must  be  less  than  0.7  in  order 
that  the  minimum  pressure  occurs  well  downstream.  Smaller  values  of  h 
are  desirable,  but,  on  the  other  hand,  we  should  like  to  make  h  as 
large  as  possible,  beca\ise  a  large  value  of  h  will  be  advantageous  in 
increasing  the  maximum  lift  coefficient  and  in  preventing  the  inception 
of  adverse  pressure  gradient  when  the  angle  of  attack  is  slightly  changed. 
Even  if  we  give  h  the  maximum  permissible  value  O.7,  the  leading-edge 
radius  amounts  to  only  60  percent  of  that  for  the  conventional  NACA  sym- 
metrical airfoil  of  the  same  thickness.   In  order  to  increase  the  leading- 
edge  radius,  it  is  required  to  increase  the  thickness,  which  in  turn  is 
accompanied  by  an  increase  in  adverse  pressiire  gradient  following  the 
minimum  pressure.  The  adverse  pressure  gradient  should  be  kept  within  a 
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certain  limit,  so  that  it  becomes  necessary  to  make  a  compromise  between 
conflicting  requirements.   Thus,  we  are  no  longer  in  a  position  to 
require  the  farthest  possible  rearward  location  of  the  minimum  pres- 
sure. We  should  also  use  a  value  of  d]_  which  is  somewhat  smaller 

than  that  mentioned  previously. 

Taking  these  requirements  into  account,  we  finally  arrive  at  the 
design  of  a  series  of  symmetrical  airfoil  sections,  the  parameters  of 
which  are  given  in  the  following  table: 


Position  of 

Section 

m 

h 

^1 

minimum  pressure 

I 

0.500 

0.35 

2.384 

0.63 

J 

.500 

.5^+ 

1.800 

.55 

K 

.1+75 

.56 

1.575 

.51 

L 

.i+50 

.58 

l.UOO 

M 

M 

.1+00 

.62 

1.150 

.37 

N 

.350 

.66 

1.000 

.21+ 

Although  section  I  is  the  most  ideal  for  delaying  the  transition,  in 
practice,  its  extraordinarily  sharp  nose  and  blunt  tail  are  drawbacks. 
On  the  other  hand,  section  N  is  too  much  compromised.   Sections  K  or  L 
seem  to  be  suitable  as  laminar-flow  airfoils  for  practical  use.   The 
ordinates  of  these  six  sections  are  given  in  table  1,  while  the  a\ixiliary 
functions  Bg  and  B^  associated  with  the  pressure  distribution  (see 

paragraph  3)  are  given  in  tables  2  and  3f   respectively.  The  shapes  of 
the  airfoil  sections  and  the  pressure  distribution  for  e  =  0.1  are 
shown  in  figure  6. 


DESIGN  OF  MEAN  CAMBER  LINE  SUTTABLE  FOR  LAMINAR-FLOW  AIRFOIIS 


8.  A  symmetrical  airfoil  set  at  zero  angle  of  attack  has  no  lift. 
In  order  to  obtain  lift,  the  center  line  of  the  symmetrical  airfoil  must 
be  curved  with  a  suitable  camber.  Since  the  effects  of  thickness  and 
camber  are  nearly  additive  with  regard  to  the  pressure  distribution,  the 
mean  camber  line  which  maintains  the  nature  of  the  pressure  distribution 
of  the  symmetrical  airfoil  will  be  such  that  it  shall  give  a  uniform  dis- 
tribution of  pressure  difference  when  the  thickness  is  removed.   Evidently, 
the  center  of  pressure  is  then  located  at  x  =  0.5,  so  that  such  a  camber 
line  has  the  drawback  that  the  travel  of  center  of  pressure  is  consider- 
able.  To  reduce  the  travel  of  center  of  pressure,  the  \miformity  of  pres- 
sure difference  should  be  satisfied  only  in  the  forward  part  of  the  chord. 
From  the  standpoint  of  designing  the  laminar-flow  airfoil,  however,  it  is 
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only  required  that  the  distribution  of  pressure  difference  is  uniform 
from  the  leading  edge  to  that  point  corresponding  to  the  mini mum  pres- 
sure of  the  symmetrical  airfoil. 

9.  When  the  angle  of  attack  a  is  small,  the  expression  for  pres- 
sure distribution  given  in  paragraph  3  may  ^^   put  into  the  form 


£=  1 


(1  +  eBg)  ±  ffAe  +  j  ^^^  aj  -  (fAg  i  eB^)a 


(' 


1  +  (fAg  1  eB^ 


f 


Since  the  effect  of  the  term  (fAg  i  eB^,!  is  very  small,  the  quantity 


fG  =  fA^ 


H 


a 


is  required  to  be  constant  in  order  that  the  camber  line  shall  not  change 
the  nature  of  the  pressiire  distribution  of  the  symmetrical  section.  The 
range  of  constancy  is  at  least  up  to  the  position  of  minimum  pressure  of 
the  symmetrical  section.   Putting  cos  I  =  u  =  2x  -  1,  and  considering 
for  simplicity  the  case  when  the  minimum  pressure  is  located  at  u  =  0, 
we  prescribe  that 


G  =  constant  =  G^ 


G  =  GqCI  -  u^)"^,  m  ^  0 


for     0  ^  X  ^  0.5,  ,  -1  ^  u  ^  0 
for     0,5  ^  X  ^  1,  0  ^  u  ^  1 


See  figure  7.  Moreover,  since 


fG 


nar 


1  -  cos  n|       1  -  cos  i 

+  a 


sin    i 


sin   I 


a  cannot  be  arbitrary,  but  must  be  so  chosen  that  the  right  side  of  the 
equation  does  not  become  infinite  at  the  leading  edge,   |  =  n.   It  is 
given  by 


8  NACA  TM  1351 


a  =  2  53  ^^ 

1 


where  (  *)  denotes  that  only  odd  integers  shoiild  he  taken  for  n.   This 

is  the  so-called  ideal  angle  of  attack  due  to  Theodorsen  (ref .  I7)  . 

Using  the  assumption  of  the  thin  wing  theory,  we  neglect  the  terms   eB^ 

and  ffAg  +  eB^j.   We  then  have 


pTx/2 

Cl  =  2(1  +  c^)fGo       a^  =    /     sin2^+l|  dl 


1  +  (1  +  m)ajn 
1+(1  +  m)(l  +  aj 


^mO  = ^L 


a  =  2  }_  nan  +  — 
1        2jt 

where  Cj^q   is  the  moment  coefficient  about  the  leading  edge  (positive 
when  nose  up)  at  C^  =  0.  Although  m  =  0  corresponds  to  the  case  of 

making  G  uniform  up  to  the  trailing  edge,  it  seems  to  be  impossible 
to  realize  a  finite  pressure  difference  at  the  trailing  edge.   Moreover, 
the  quantity  -Cmo/CL  (which  represents  the  degree  of  center  of  pres- 
sure travel)  is  as  large  as  0.25  in  this  case.   If  m  >  0,  G  vanishes 
at  the  trailing  edge,  and  -CmO/^L  decreases  as  m  increases,  tending 
to  0  as  m  approaches  c».   m  =  00  corresponds  to  the  case  when  G  =  0 
in  the  rear  half  of  the  chord.   Increasing  the  value  of  m,  however, 
steepens  the  pressure  gradient,  so  the  value  of  m  from  3  "to  5  seems  to 
be  adequate. 
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Now,  since 


fG  sin  1  =  2  2 na^  cos  n^  -  [  2  2 na^  +  —  j  cos  I  +  — 

1  \   1        2n/        2n 


the  slope  of  the  camber  line  having  the  prescribed  distribution  of  G 
is  given  by 


dM 


dx  1  sm 


sin  nl  1     P'^        r^  _  dl' 

;;—  =   -  ~   /       2   2 na^  cos  nl     j- 

rt  Jo  1  cos    I   -  cos    I 


2jt 


1   + 


1  +  a 


m 


-,  u  /,  o\m  1  -  u 

log  +   (1   -   u^)      log  

1  +  u 


u 


1    (1   -   u2)"   -    (1   -   v2) 


m 


dv 


U    -    V 


+  2   y>      na. 


"n 


The  ordinate  of  the  camber  line  is  obtained  by  the  integration 


M  =  f  ^  dx 
Jo  <ix 


2 na^^  may  be  determined  by  the  condition  that  M  =  0  at  x  =  1.   We 

1 
call  Dju  the  camber  line  thus  determined.  The  equations  for  camber 

lines  for  m  =  0,  1,    3,   5^  and  <»,  namely  Do,  Di,  D3,  D^,   and  D„, 


10 
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are  given  below,  their  important  characteristics  being  simmarized  in  the 
table.   It  is  to  be  noted  that  f  is  the  maximum  value  of  M,  e   is 
the  absolute  value  of  the  zero-lift  angle,  and  e   and  a  are  measured 
in  radians. 


m 

cx/Cl 

^/Cl 

f/CL 

-CmO/CL 

0 

0 

0.1592 

0.0552 

0.2500 

1 

.0380 

.1211 

.0711 

.1750 

^ 

.0609 

.0983 

.0790 

.1213 

5 

.0703 

.0888 

.0816 

.0979 

00 

.1103 

.Oi^88 

.087^+ 

0 

M 


=  1 


2  log 


^^—  |(1  -  u)  log  (1  -  u)  +  (1  +  u)  log  (1  +  u)  ^ 
og  2  [_  J 


D. :  Jl  M  =  i  ( 5  +  u)  log  2  -  -  ( 1  +  u)  log  ( 1  +  u)  +  i-  u3  log  I  u 
Cl     5  5  5 


i  (1  -  u)2(2  +  u)  log  (1  -  u)  +  i  (1  -  u2) 
5  5 


Dn:   _il  M  =  —  (51  +  19u)  log  2  -  —  (1  +  u)  log  (l  +  u)  + 
(-L     51  51 

—  u3(35  -  21u2   +  5u^)  log    lu  I  - 
51        \  '  ' 

-i-  (1  -  u)^(l6  +  29u  +  20u2  +  5u3)    log  (1  -  u)    + 
51 


-^  (1  -  u2)(i76  .-  8lu  -  172u2  +  30u3  +  60u^) 
612 
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D.^:      il  M  =  J^  (9if9  +  l+37u)    log  2   -  221  (1  +  u)    log   (l  +  u)    + 
?       Cl  9^9  9^9 

-—-  u3 (1155  -  i386u^  +  990u^  -  385u^  +  63u^)  log  | u  |  - 

-^  (1  -  u)^(256  +  8U3U  +  I2l8u2  +  938u3  +  378u^  + 

63u5)  log  (1  -  u)   4  — - —  (1  -  u2) (35072  -  28535U  -  66o88u2 

113880  ^ 

3l680u3  +  68792u^  -  17J4-30u5  -  36l20u^  +  3780u'i'  +  7%0u°) 

\'   T  =  -^—    |(1  +  ^)    l°g  2   -    (1   +  u)    log   (1  +  u)    +  u  log   |u|  [ 
1       log  3     I  -> 


Shapes  of  these  camber  lines  are  shown  in  figure  7,  "the  ordlnates 
of  them  are  given  in  table  k,   and  the  auxiliary  functions  A^  and  Ag 
(see  paragraph  3)  and  the  pressure  difference  distribution  G  are  given 
in  tables  5,  6,  and  J,   respectively. 

10.  The  calculation  made  previously  is  only  approximate,  neglecting 
the  thickness.   It  is  therefore  desirable  to  check  the  result  by  actually 
calculating  the  pressure  distribution  for  the  specified  angle  of  attack 
taking  both  camber  and  thickness  into  account.   As  an  example,  we  con- 
struct an  airfoil  by  applying  the  thickness  form  K  with  e  =  O.I5  nor- 
mal to  the  chord  around  the  camber  line  D^  with  f  =  0.02   (the 
resulting  airfoil  is  designated  as  DcK  -  2OI5).   We  calculate  the  pres- 
sure distribution  by  the  formula  of  paragraph  3  for  the  optimum  design 
condition  a  =  0.99°  and  C^  =  0.2il5.  The  result  is  shown  in  figure  8. 

The  nature  of  the  pressure  distribution  remains  similar  to  that  of  the 
symmetrical  airfoil,  so  we  may  consider  that  the  approximate  determina- 
tion neglecting  thickness  gives  results  sufficiently  accurate  for  practi- 
cal purposes. 

11.  In  designing  the  camber  line     B^,   we  have  assiraied  for   simplicity 

that  the  pressure  difference  G   is  constant  for  u  "S  0.  This  corresponds 
to  the  case  when  the  symmetrical  airfoil  has  its  minimum  pressure  in  the 
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neighborhood  of  u  =  0.  Therefore,  the  camber  line  Dm  is  adequate  to 
be  combined  with  the  symmetrical  section  J  or  K.  If,  however,  the 
symmetrical  section  is  adopted  in  which  the  minimum  pressure  is  located 
further  upstream,  it  is  not  only  not  necessary  to  maintain  G  constant 
up  to  u  =  0,  but  also  of  disadvantage  because  it  makes  it  difficult  to 
reduce  the  value  of  -Cjuq/CL' 


To  reduce  the  range  over  which 
we  may  proceed  in  the  following  way. 

should  be  constant  from  u  =  -1  to 

u-i  =  —  (1  +  3^),   we  prescribe  that 
4 


G   should  be  maintained  constant. 
Assuming  for  Instance  that  G 

u  =   -  -7,   and  xising  a  new  variable 


G  =  Gr 


for 


-1^  u^   - 


1 
3' 


-|^ui<=0 


G   =  Go(l   -   ui2) 


2^3 


for 


is  u^  1, 


0  ^  u^  ^  1 


The  calculation  may  be  performed  similarly  to  the  case  of  D 


m* 


The 


resulting  camber  line  is  designated  as  Fo.  The  camber  line  lying  in 
the  middle  between  Fo  and  Dc  is  also  designed,  and  designated  as  Ej^. 

Their  important  characteristics  are  given  in  the  following  table  together 
with  those  of  D^.   Other  niimerical  data  for  these  camber  lines  are  given 

in  tables  k   to  7« 


Camber  line 

a/CL 

€/Cl 

f/CL 

-•^mO/CL 

^5 
^3 

0.0703 
.0752 

.0764 

0.0888 
.081+0 
.0827 

0.0816 

.0813 
.0795 

0,0979 
.0859 
.0833 

Since  the  camber  lines  Ei^  and  Fo   enable  us  to  maintain  G  con- 
stant up  to  the  point  x  =  0.i<-2  and  x  =  O.33,  respectively,  they  are 
adequate  to  be  combined  with  the  symmetrical  sections  L  and  M,  respec- 
tively. The  pressvire  distribution  is  shown  in  figure  10  for  the  airfoil 
obtained  by  applying  the  thickness  form  M  with  e  =  O.I5  around  the 
camber  line  F-^  with  f  =  0.02.  The  optimum  design  condition  corresponds 

to  a  =  1.10°  and  C^  =  0.252. 
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EXPERIMENTS  ON  LAMINAR -FLOW  AIKFOII£ 


12.  In  order  to  ascertain  whether  it  is  possible  to  prevent  the 
forward  movement  of  the  boimdary  layer  transition  "by  shifting  the  mini- 
mum pressure  on  the  airfoil  surface,  we  have  to  perform  experiments  in 
a  low  turbulence  wind  tunnel  or  on  the  actual  airplane  in  flight.  When 
the  Reynolds  number  is  not  too  large,  however,  we  can  still  use  a  con- 
ventional wind  tionnel  in  which  the  stream  turbulence  is  relatively  small. 
So  we  made  at  first  comparative  measurements  on  two  symmetrical  airfoils, 
NACA  0010  and  L.B.  2k   in  the  1.5  m  wind  tunnel  of  the  Aeronautical 
Research  Institute.   L.B.  2k   is  a  laminar-flow  airfoil  of  10  percent 
thickness,  already  shown  in  figure  3«  The  theoretical  pressure  distri- 
bution is  also  given  in  figure  11.  The  minimum  pressure  is  located  at 

X  =  0.6k,   and  the  laminar  separation  at  x  =  0.77.  The  wind  timnel  was 
of  the  lowest  turbulence  level  available  for  the  author,  the  critical 

Reynolds  number  of  a  sphere  being  3 '66  X  lO-'  and  the  transition  Reynolds 

number  of  a  flat  plate  I.05  X  lo3  (see  paragraph  28) .   In  order  to  raise 
the  Reynolds  number  as  high  as  possible,  unusually  large  models  were 
used.  They  were  made  of  laminated  mahogany,  of  highly  polished  surface, 
of  0.8  m  span,  of  1.2  m  chord,  and  fitted  with  end  plates  1.3  m  X  0.6  m. 
Since  the  model  was  large  compared  with  the  size  of  the  tunnel  and  the 
end  plates  were  not  sioff iciently  large,  the  results  for  a  given  airfoil 
may  not  correspond  even  approximately  with  those  for  the  same  airfoil  in 
an  undisturbed  two-dimensional  flow.   Our  object,  however,  was  merely 
to  ascertain  the  relation  between  pressure  distribution  and  transition, 
and  it  seemed  reasonable  to  expect  that  the  relation  will  not  be  seriously 
affected  by  limitations  in  the  conditions  of  the  experiments.  As  a  matter 
of  fact,  marked  difference  was  found  in  the  calculated  and  measured  dis- 
tributions of  pressure,  the  latter  of  which  was  measured  along  the  median 
section  of  the  model  with  a  static  tube  of  1  ram  diameter'  (fig.  12). 2 
This  discrepancy,  however,  is  immaterial,  since  our  object  was  merely  to 
compare  the  two  airfoils,  both  of  which  are  affected  quite  similarly  by 
experimental  limitations. 

13.  The  angle  of  attack  of  the  model  was  zero,  and  the  wind  speed 
was  varied  from  6  to  kO   m/s.  The  local  drag  of  the  median  section  was 
determined  from  wake  measurements,  that  were  made  in  the  section  11  cm 
behind  the  trailing  edge.  Measurements  of  static  and  total  press\ires  in 
the  wake  were  made,  respectively,  with  a  static  tube  of  2.5  ram  external 
diameter  and  a  pitot  tube  with  a  flattened  mouth  of  O.65  mra  external 
depth  and  2.6  ram  width.  The  profile  drag  coefficient  Cdq  was  obtained 


P  6 

The  measured  values  are  those  for  a  Reynolds  number  of  about  2  X  10°. 

The  distribution  of  pressure  changes  but  little  with  the  Reynolds  number. 
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from  the  measured  pressiires  tj   Jones'  formula  (ref.  I8) .  Figure  I3  3 
presents  C^q  plotted  against  Reynolds  number  R  referred  to  chord 

length.  For  a  lower  range  of  R,  the  drag  of  L.B.  2k   is  higher  than 
that  of  NACA  0010,  the  reason  probahly  being  that  a  turbvilent  boundary 
layer  associated  with  the  laminar  separation  is  established  at  a  higher 
Reynolds  number  for  the  former  airfoil  than  for  the  latter.  For  a  higher 
range  of  R,  however,  the  condition  is  reversed,  L.B.  2U  giving  a  drag 

less  than  half  that  of  NACA  0010  for  R  higher  than  2  X  10^.  This  is 
probably  due  to  the  fact  that  the  transition  may  occur  much  later  for 
L.B.  2k   than  it  does  for  NACA  0010,  as  also  observed  from  the  compari- 
son of  wake  conditions  for  the  two  airfoils  (fig.  1^). 

ik.    In  order  to  verify  the  aforementioned  supposition,  a  pitot  tube 
with  a  flattened  mouth  of  external  depth  O.9  mm  and  width  2.7  mm  was 
placed  in  contact  with  the  airfoil  surface,  and  the  wind  speed,  and  con- 
sequently the  Reynolds  n\amber  R,  were  determined  at  which  the  indicated 
total  pressure  G*  divided  by  the  dynamical  pressure  q  of  the  undis- 
turbed stream  begins  to  rise  suddenly.  The  results  are  shown  in  fig- 
ure 15.  From  this  figure,  the  dependence  of  the  transition  point  on 
Reynolds  number  as  shown  in  figure  I6  is  obtained.  At  the  same  Reynolds 
niomber,  the  transition  occurs  much  farther  from  the  leading  edge  for 
L.B.  2k   than  for  NACA  0010.  Even  at  the  highest  Reynolds  number  reached, 
L.B.  2k   has  a  transition  as  far  back  as  x  =  O.8O.   This  is  somewhat 
beyond  the  laminar  separation  point,   x  =  0.11,   which  is  calculated  from 
the  theoretical  pressure  distribution.   However,  this  is  not  contradictory, 
because  the  actual  pressure  distribution  differs  from  the  theoretical  one 
in  a  manner  to  delay  the  transition  (fig.  12). 

15.  With  further  increase  in  Reynolds  number,  the  transition  may 
move  toward  the  leading  edge,  but  it  seems  improbable  that  the  transition 
moves  forward  beyond  the  minimum  pressure.   It  is  highly  desirable  to 
check  this  point  also  by  wind  tunnel  experiments,  but  all  the  wind  tunnels 
now  available  to  the  author  are  of  no  use  for  making  measurements  at 
sufficiently  high  Reynolds  numbers,  becaiise  the  transition  is  prematurely 


3ln  this  figure,  the-  ciirves  L  and  T  represent  the  drag  of  a 
flat  plate  when  the  boundary  layer  is  entirely  laminar  and  entirely 
turbulent,  respectively.   The  curves  NV  and  irF  represent  the  drag 
of  airfoil  NACA  OOO9  measured  in  the  NACA  Variable-Density  Wind  Tunnel 
and  NACA  Full-Scale  Wind  Tunnel,  respectively. 

In  this  figure  G  and  p  are  the  total  and  static  press\ires  in 
the  wake,  respectively,  and  Gq   is  the  total  pressure  outside  the  wake, 
all  being  measured  from  the  static  pressure  of  the  undisturbed  stream, 
y  is  the  distance  across  the  wake,  and  t   is  the  chord  of  the  model. 
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induced  "by  the  turbulence  of  the  stream  (see  paragraphs  27  to  29) •   It 
seems  urgent  to  build  a  special  wind  tunnel  of  low  turbulence  level. 
For  the  present,  however,  it  is  simplest  to  rely  upon  experiments  in 
actual  flight.  Such  a  hope  of  the  author  was  fortunately  realized  by 
the  specially  planned  flight  experiment,  which  was  performed  at  the 
Navy  Aeronautical  Technical  Arsenal  (ref.  19). 

16.  The  airplane  used  for  the  experiment  was  a  biplane;  two  portions 
of  the  lower  wing,  each  of  1.1  m  span,  were  covered  with  the  airfoil  to 

be  tested.   The  test  portions  were  of  chord  2.^  m,  made  of  Japanese  Hinoki, 
highly  polished,  and  fitted  with  a  partition  fence  of  small' height  at 
both  ends.  Two  test  portions  were  placed  symmetrically,  pressure  distri- 
bution and  wake  measurements  being  performed  on  the  starboard  portion, 
while  the  boundary  layer  was  observed  on  the  port  portion.   The  airfoil 
section  was  not  one  of  the  most  appropriate  design  now  considered,  because 
it  was  required  to  put  it  on  the  original  section  of  the  airplane,  and, 
moreover,  to  determine  the  section  before  completion  of  the  final  design 
calculation.   It  has  the  following  characteristics: 

Mean  camber  line:   M  =  0.0667x(l  -  x)(l  -  x  +  x^),  f  =  0.0125 

Thickness  distribution:   e  =  0.12,  m  =  0.1+5,  h  =  O.56,  6.^   =  I.60 

The  camber  line  is  similar  to  Dq  of  paragraph  9^  tut  there  exists  a 

slight  lack  of  uniformity  of  G  in  the  neighborhood  of  leading  and 
trailing  edges.  The  thickness  distribution  is  similar  to  L  of  para- 
graph 7>  "but  the  trailing  edge  slope  is  somewhat  larger  than  L. 

17.  Results  of  flight  experiments  are  summarized  in  figure  17 
and  18.   In  figure  I7,  the  section  lift  coefficient  Cl,  obtained  by 
integrating  the  pressure  distribution  curve,  is  shown  by. a  broken  line 
plotted  against  the  Reynolds  number  R  referred  to  the  flight  speed 

and  chord  length,  and  C^  is  again  shown  by  a  solid  line  plotted  against 

the  profile  drag  coefficient  Cj)_  determined  from  the  wake  measurements. 

In  figure  I8,  the  measured  pressure  distribution  is  shown  in  comparison 
with  the  theoretical  one  (two-dimensional  potential  flow)  having  the 
same  value  of  Cl.  The  transition  points  estimated  from  the  change  in 
boundary  layer  velocity  profiles  are  also  marked.  Generally  speaking, 
the  measured  pressure  distribution  agrees  fairly  well  with  the  theo- 
retical one,  although  a  slight  difference  appears  when  C^  becomes 
large.  An  adverse  pressure  gradient  is  found  on  the  lower  surface  when 
Cl  is  small,  thus   resulting  in  the  transition  point  being  observed 
^mexpectedly  far  forward.   Such  a  discrepancy  in  pressure  distribution 
as  compared  with  the  theoretical  one  seems  to  be  probably  due  to  the 
fact  that  the  span  of  the  test  portion  was  not  sufficiently  large.  As 
a  result,  the  profile  drag  coefficient  C-q       has  the  minimum  value  0.00^2 
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at  about  Cl  =  O.26,   which  is  larger  than  the  value  Cl  =  0,l8  theo- 
retically estimated  on  the  assumption  that  the  transition  occurs  far 
hack  on  hoth  upper  and  lower  surfaces.   Therefore,  the  observed  value 
°f  '^DO'  al'though  much  smaller  than  that  of  the  conventional  airfoils, 

seems  to  be  still  somewhat  large  when  compared  with  the  optimum  case. 
At  any  rate,  however,  no  transition  was  found  to  occur  upstream  of  the 
minimum  pressure.   It  is  important  to  note  that  such  experimental  evi- 
dence was  obtained  on  an  airfoil  section  in  which  the  minimum  press\are 
is  located  further  downstream  than  on  the  conventional  one.  This  finding 
will  give  valuable  data  to  establish  a  basis  for  design  of  the  laminar- 
flow  airfoils. 


ESTIMATION  OF  THE  DRAG  OF  lAMINAR-FLOW  AIRFOIIS 


18.  As  mentioned  previously,  the  results  of  flight  experiments 
seem  to  support  the  basis  for  the  design  of  laminar-flow  airfoils, 
namely,  the  possibility  of  maintaining  the  boimdary  layer  laminar  at 
least  up  to  the  minimum  pressure  point.   It  is  interesting,  therefore, 
to  estimate  the  drag  of  laminar-flow  airfoils  by  assuming  a  laminar 
boundary  layer  from  the  leading  edge  to  the  minim-urn  pressure  point  and 
a  turbulent  boundary  layer  downstream  to  the  trailing  edge. 

For  the  laminar  boundary  layer,  the  momentum  thickness  is  given  by 


2   0.l44v 


Ui 


ds 


with  a  sufficient  approximation  (ref.  16)  ,  where  u-^      is  the  velocity 

outside  the  boimdary  layer  and  s   is  the  distance  measured  along  the 
airfoil  surface  from  the  forward  stagnation  point.   Writing  t  for  the 
chord  length  and  V  for  the  velocity  of  the  undisturbed  flow  (velocity 
of  flight) ,  and  putting 


u-L  =  UV 


s  =  at 


R=Xt 


we  have  the  nondimensional  expression 
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a 
u2 


O.kk 
RU„ 


(^)' 


da 


where  the  subscript  a  refers  to  the  point  of  minimum  pressure. 
Applying  then  the  solution  due  to  Buri  (ref .  20)  for  a  turbulent 
boundary  layer  assumed  to  extend  from  the  minimum  pressure  point  to 
the  trailing  edge,  we  have  the  result 


'e\W\ii/h 


^) 


u 


=  0.0162R 


■lA  r^  ^ 


da 


where  the  subscript  b  refers  to  the  trailing  edge.  The  niunerical 
values  originally  given  by  Buri  are  slightly  modified  so  as  to  agree 
with  measurements  when  applied  to  the  flat  plate. 

According  to  Squire  and  Young  (ref.  2l) ,   the  profile  drag  coeffi- 
cient is  given  by 


'Do 


2 
t 


(^u 


az)u. 


3.2 


where  the  subscripts  u  and  Z   refer  to  the  upper  and  lower  surfaces, 
respectively.   The  exponent  3.2  of  U-j^  has  been  obtained  by  assuming 

the  ratio  of  displacement  and  momentum  thicknesses  equal  to  l.U.  But  the 
ratio  seems  to  exceed  l.h   near  the  trailing  edge,  so  we  replace  3'2  by 
3.^  with  a  view  to  improving  the  accuracy  and  at  the  same  time  to 
simplifying  the  algebra.   Since 


T^Ub 


Z.h 


ea\5A„  itA 


(^) 


u= 


+  0.0162R 


-lA 


'b  ,,k 


U     da 


V5 
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we  have 


Cdq  =  0.07i+R 


-1/5 


(l,R-3/8  .  T>/5  .  (l,H-3/8  ,  ,^)V5 


where 


L  =  3m 


9/8 


a  I, 
U^  da 


5/8 


T  = 


^«^b 


lA  da 


If  the  velocity  distribution  U3_  =  UV   is  calculated  by  assuming  the 

potential  flow  of  an  ideal  fluid,  it  is  desirable  to  modify  the  distri- 
bution to  take  accovmt  of  the  effect  of  separation  near  the  trailing 
edge.  We  tentatively  modified  the  distribution  of  U  such  that  the 

Buri  parameter  T   =  (0/U]^)(dui/dsWu2.9/v)  '    at  the  trailing  edge  for 

the  case  when  the  boundary  layer  is  assumed  turbulent  from  the  leading 
edge,  namely 


0.0081  /dU'^ 


U> 


Vda 


U^  da 


shall  not  become  smaller  than  -O.O6.   In  almost  all  the  cases,  values 
of  U-^  thus  modified  are  found  in  the  range  between  0.95  and  1.00. 

19.  Applying  this  method  of  calculation,  the  profile  drag  coeffi- 
cient Cj)   is  estimated  first  for  a  series  of  symmetrical  airfoils  set 

at  zero  angle  of  attack.   The  series  consists  of  the  six  symmetrical 
airfoils,  I,  J,  K,  L,  M,  N,  as  given  in  paragraph  7  and  the  NACA  con- 
ventional airfoil.   Values  of  Cp   at  R  =  2  X  lo'''  for  three  different 

thicknesses  (maximum  thickness  in  terms  of  chord  e  =  0.10,  O.I5,  0.20) 
are  shown  in  figure  I9  plotted  against  the  position  of  minimum  pressiire. 
Cj)q  seems  to  decrease  almost  linearly  as  the  minimum  pressure  is  shifted 

backward,  the  most  ideal  airfoil  I  giving  a  value  about  half  of  that  of 
the  NACA  conventional  airfoil.   If  it  is  desired  to  realize  a  profile 
drag  of  two-thirds  of  the  conventional  airfoil,  it  will  be  required  to 
use  the  symmetrical  airfoil  L  with  the  maximum  thickness  located  at 
h^   percent  chord  from  the  leading  edge. 
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Effect  of  camber  is  relatively  small.   If,  for  instance,  the 
center  line  of  the  symmetrical  airfoil  K  with  e  =  O.I5  is  curved 
into  the  camber  line  D5  with  f  =  0.02   (see  paragraph  9)  >  "the  esti- 
mated increase  in  profile  drag  at  the  optimum  angle  of  attack  is  only 
0.0001.  For  f  =  O.Oij-,  it  is  O.OOO3. 

Finally,  we  compare  the  laminar-flow  airfoil  with  the  most  exten- 
sively used  airfoil,  NACA  23012,  for  which  the  leading-edge  radius  is 
0.0158,  and  the  optimum  lift  coefficient  corresponding  to  the  minimvun 
profile  drag  coefficient  is  about  O.I5.   If  we  consider  the ' symmetrical 
airfoil  section  K   (h  =  0 . 56)  combined  with  the  camber  line  Dr,  it  is 

necessary  to  use  the  thickness   e  =  O.15  in  order  to  obtain  the  same 
magnitude  of  leading-edge  radiiis,  and  the  camber  f  =  0.012   in  order  to 
realize  the  optimum  lift  coefficient  0.15.^  Therefore  we  construct  an 
airfoil  by  applying  the  symmetrical  form  K  with  e  =  O.I5  normal  to 
the  chord  around  the  camber  line  D5  with  f  =  0.012.  We  call  it 
D5K  -  1215.  The  angle  of  attack  corresponding  to  C^  =  0.15  is  1.56° 
for  NACA  23012  and  0.60°  for  D5K  -  1215-   The  pressure  distribution  for 
that  condition  is  shown  in  fig\ire  20. 

We  then  estimate  C-Qq     for  the  two  airfoils  by  the  method  explained 

previously.  The  results  are  shown  by  broken  lines  in  figure  21.   In 
order  to  check  the  results,  measured  values  taken  from  various  sources 
for  the  two  airfoils  and  similar  airfoils  are  also  plotted  in  the  same 
figure  by  different  marks.   The  mark  o  refers  to  the  value  obtained 
by  flight  experiments  on  a  smooth  surface,  and  •  refers  to  that 
obtained  by  wind  tunnel  experiments  where  the  stream  turbulence  has  no 
effect  on  transition.  The  mark  +  refers  to  the  flight  experiment  on 
a  rough  surface,  while  X  refers  to  the  wind  tunnel  experiment  where 
the  stream  turbulence  causes  the  transition  to  occur  prematurely.   There- 
fore, only  o  and  •  are  adequate  for  our  present  purpose.   Drawing 
ciirves  through  these  points  and  extrapolating  to  higher  Reynolds  numbers, 
we  find  that  the  result  agrees  fairly  well  with  the  estimated  values. 
Therefore,  we  may  consider  that  the  method  of  estimating  Cj)_   is  suffi- 
ciently accurate  at  the  Reynolds  numbers  corresponding  to  actual  flight 

5The  calculation  developed  in  paragraphs  8  to  11  refers  to  the  poten- 
tial flow  of  an  ideal  fluid,  so  that  it  gives  the  slope  of  lift  curve 

dCr 

=  2n.   In  real  fluids,  however,  the  slope  of  lift  curve  amounts  to 

da 

only  80  to  90  percent  of  the  theoretical  value.   If  we  take  this  effect 

into  accoinit,  we  have  to  increase  the  necessary  amount  of  f  by  10  to 
20  percent  in  order  to  realize  the  given  lift  coefficient.   However, 
such  a  slight  change  in  the  value  of  f  will  scarcely  affect  the  esti- 
mation of  Cj)„. 
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conditions.  Comparison  of  two  airfoils,  laminar-flow  and  conventional, 
also  suggests  the  possibility  of  kO   percent  reduction  in  profile  drag 
by  using  a  fairly  practical  laminar-flow  airfoil. 


CONSIDERATION  OF  THE  AIRFOIL  WITH  UNIFORM  DISTRIBUTION  OF  PRESSURE 


20.  The  fact  that,  so  far  as  flight  experiments  with  smooth  wings 
are  concerned,  the  boundary  layer  transition  occurs  only  in  the  region 
of  rising  pressure,  not  only  warrants  the  principle  of  designing  the 
laminar-flow  airfoil  by  shifting  the  minimum  pressure  backward,  but 
also  suggests  the  possibility  of  delaying  the  transition  by  using  an 
airfoil  with  uniform  distribution  of  pressure.  Therefore,  in  para- 
graphs 21  to  23,  the  shape  of  such  a  symmetrical  airfoil  is  determined 

by  a  method  similar  to  that  used  for  designing  the  camber  line  of  laminar- 
flow  airfoils,  and  the  airfoil  was  examined  by  wind  tunnel  experiments. 
In  paragraphs  2k   to  26,  a  calculation  is  made  to  inquire  about  the  method 
of  sucking  away  the  boundary  layer  over  the  region  of  rising  pressxire  in 
such  a  way  that  the  boundary  layer  velocity  profile  shall  remain  the  same 
as  that  for  the  point  of  minimum  pressure. 

21.  Consider  the  syrmnetrical  airfoil  set  at  zero  angle  of  attack. 
According  to  the  formula  of  paragraph  3^  "the  pressure  distribution  is 
given  by 


€•=  1 


sin  nl 


1  +  2  X;^  nb^ 

1      sin  i 


1  + 


CO  t 

T COB  n| 

sUnb^  - 

1     sin  I 


where  p  is  the  pressiire  acting  on  the  airfoil  surface,  measured  from 
the  static  pressiire  of  the  undisturbed  stream,   q  is  the  dynamic  pres- 

Bure   of  the  imdisturbed  stream,  x  =  ^  (l  +  cos  i)      is  the  coordinate 

along  the  chord,  and  the  half -thickness  of  the  airfoil  is  expressed  in 
the  form 


T  =  ^  bjj  sin  ni 
1 
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If  the  thickness  is  sufficiently  small,  the  square  of 


dx 


-2  Yi   nbn 

1 


cos  nl 


sm 


i 


may  be  neglected,  so  that  the  condition  of  uniform  distribution  of  pres- 
sure is  satisfied  by  putting  all  the  coefficients  b^,  other  than  bj^, 

equal  to  zero,  namely,  by  an  elliptic  section.   In  order  to  take  the 
thickness  into  account  approximately,  we  substitute  the  value  of  dT/dx 
for  the  elliptic  section  into  the  denominator  of  the  expression  for  p/q. 
Then,  writing  e  for  the  maximum  thickness  in  terms  of  the  chord,  we 
get 

00 

sin  I  +  2  2_  nbjj  sin  n| 

"^  .  2t    2   2t 

sm  5  +  e  cos  5 


hence 


2  2 nbjj  sin  n4  =  -  sin 


i  +   BVsin^^  +  e^cos^l 


1/2 
where  B  is  the  constant  value  of  (l  -  p/q)    .  We  have  therefore 
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^i'^  ^   S  =   -  T    /      2   ^  nb^  sin  nl 


dx  K 


^0  1 


cos  i   -  cos  ^' 


jt 


'0 


>, 


in   r    -  B 


sin 


Jsin2^'    +  e^cos^l' 


Bin  I'    dl' 
cos  i   -  cos  ^' 


^ 


B  ^ 
u  -  —  •< 


2ku  cos     e   + 


ll   -  k^u^   log 


(1  +  u)  (1   -  k^u  +  eflT^V)    I 
(1  -  u)(l   +  k^u  +  eU   -  k^u^) 


where 


u  =  cos  I    =  2x  -  1  k  =  \|: 


1  -  e^ 


Upon  integrating  we  get 


T  =  <||k  cos-^e  -  iHl  -  u2 


_B 
2 


u 


^ii\| 


1  -  k^u^     (1  +  u)  (1  -  k^u  +  e^l  -  k^u^) 
log 


du 


1  -  u^ 


(1  -  u)(l  +  k^u  +  e\|l  -  k^u^j 


The  integral  is  evaluated  by  a  numerical  method,  and  the  value  of  the 
constant  B  determined  from  the  condition  that  T  =  ie  when  u  =  0. 

The  numerical  results  for  three  values  of  e  are  given  in  table  8,  T 
and  dT/dx  being  expressed  in  terms  of  those  for  the  elliptic  section. 
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The  shapes  of  the  airfoils  are  shown  in  figure  22.  The  shape  resembles 

an  elliptic  section,  although  it  is  somewhat  fuller  at  the  ends.   It  will 

also  he  seen  that  the  constant  B,  as  shown  in  the  following  table,  is 

1/2 
slightly  smaller  than  1  +  e,  the  maximum  value  of  (l  -  p/q)    for  the 

elliptic  section.   Since  the  values  of  (dT/dx)  ^   (-eu/l|l  -  u^ )  are  not  far 
different  from  1,    it  seems  to  be  sufficiently  accurate  to  substitute 
the  value  of  dT/dx  for  the  elliptic  section  into  the  denominator 
of   1  -  p/q. 


e 

B 

0.1 
.2 
.3 

1.097 
1.188 
1.273 

22.  The  uniform  distribution  of  pressure  requires,  however,  an 
infinite  pressure  gradient  at  both  leading  and  trailing  edges.   In  order 
to  see  to  what  degree  such  a  sharp  pressure  gradient  may  be  realized  in 
actual  fluids,  measurements  were  made  on  a  model  of  the  airfoil  section 
with  uniform  distribution  of  pressure  with  e  =  0.1  (we  call  it  U.P.  OOIO) 
in  the  1.5  m  wind  tunnel  of  the  Aeronautical  Research  Institute.  The 
model  was  made  of  laminated  mahogany,  of  0.8  m  span,  of  0.8  m  chord,  and 
fitted  with  end  plates  1.3  ni  X  0.6  m.  Measurements  of  pressure  distribu- 
tion, wake  traverse  and  boiondary  layer  transition  were  similar  to  those 
already  mentioned  in  paragraphs  12  to  lU. 

The  pressure  distribution  along  the  chord  is  shown  in  figure  23 
for  three  values  of  R,  the  Reynolds  number  referred  to  chord  length. 
The  observed  value  is  somewhat  high  compared  to  the  theoretical 

value  —  =  -O.2O3,  the  discrepancy  probably  being  due  to  the  excessive 

size  of  the  model  in  proportion  to  that  of  the  wind  tunnel.  At  any  rate, 
however,  the  pressure  distribution  is  nearly  uniform.   The  lack  of  uni-  , 
formity  exists  at  both  edges  due  to  the  impossibility  of  realizing  the 
infinite  pressure  gradient.   The  boundary  layer  separates  near  the 
trailing  edge,  but  the  effect  of  separation  becomes  small  as  the  Reynolds 
number  increases.  This  scale  effect  seems  to  be  of  the  same  nature  as 
that  responsible  for  the  sudden  drop  in  sphere  drag;  the  boundary  layer 
separates  in  a  laminar  state  when  the  Reynolds  number  is  low,  while  it 
becomes  turbulent  before  separation  when  the  Reynolds  number  is  high, 
thus  being  able  to  proceed  against  a  larger  pressure  gradient.   This  is 
also  seen  from  the  measurements  in  the  wake,  where  the  indentation  of  the 
curve  of  total  pressure  distribution  is  shallow  and  wide  for  low  Reynolds 
numbers,  while  it  becomes  deep  and  narrow  as  the  Reynolds  number  increases. 
As  a  result,  the  profile  drag  coefficient  C^q  decreases  considerably 
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as  the  Reynolds  number  increases,  as  shown  in  figure  2k.      Transition  to 
turbulence  was  found  very  near  to  the  trailing  edge,  occurring  downstream 
of  X  =  0.9  in  the  range  of  measurements.   Measured  values  of  the  drag 
of  the  model  when  a  piano  wire  of  0.5  inm  diameter  was  placed  at  x  =  0.8 
and  X  =  0.9,  respectively,  are  also  plotted  in  the  same  figure.  The 
drop  in  drag  occurs  at  a  lower  value  of  the  Reynolds  number  when  the 
surface  is  roughened  by  the  wire. 

23.  The  profile  drag  coefficient  of  the  airfoil  U.P.  0010  is  shown 
below  in  comparison  with  other  symmetrical  airfoils  at  R  =  2.2  x  10°: 

NACA  0010  Cj)  =  0.006i+  (fig.  I3) 

L.B.  2k  "   0.0032  (fig.  13) 

U.P.  0010  0.0059 

U.P.  0010,  a  wire  at  x  =  O.9  0.004J+ 

All  the  airfoils  are  of  10  percent  thickness.  NACA  0010  is  a  conventional 
airfoil,  and  L.B.  2k   is  a  laminar-flow  airfoil  with  far  back  minimum 
pressure.  The  drag  of  U.P.  0010  is  between  that  of  these  two  airfoils, 
the  drag  when  a  wire  is  placed  being  nearly  the  mean  of  the  two.  This 
result  seems  to  be  interesting  in  that  the  drag  of  an  airfoil  with  a 
blimt  tail  is  smaller  than  commonly  considered. 

The  airfoil  with  uniform  distribution  of  pressure  will  also  probably 
be  favorable  when  used  at  high  subsonic  speeds.   Even  if  the  shock  wave 
occurs  at  high  subsonic  speeds,  the  increase  in  drag  will  remain  small 
when  the  boundary  layer  does  not  separate.   This  expectation  was  really 
verified  by  the  experiment  due  to  Kawada  and  Kawamiira  (ref.  22),  the 
drag  of  the  airfoil  U.P.  0010  being  smaller  at  high  Mach  numbers  as 
compared  with  other  airfoils. 

2k.   From  the  fact  that  the  boundary  layer  transition  occurs  only  in 
the  region  of  rising  pressure,  we  may  also  expect  suction  of  the  boundary 
layer  to  delay  transition.  For  example,  if  the  boundary  layer  is  sucked 
into  a  slot,  there  is  a  well-known  sink  effect  (ref.  23)  which  relieves 
the  adverse  pressure  gradient  somewhat  upstream  of  the  slot.  We  may  con- 
sider an  alternative  possibility.  That  is,  we  assume  that  the  boundary 
layer  is  laminar  in  the  region  of  falling  pressure,  and  that  it  remains 
laminar  also  in  the  region  of  rising  pressure  provided  that  the  boimdary 
layer  profile  is  the  same  as  that  at  the  minimum  pressure  point.   We 
then  ask  what  suction  arrangement  must  be  applied  in  order  to  realize 
such  a  condition. 

25.  We  denote  by  s  the  coordinate  measured  along  the  surface, 
y  perpendicular  to  the  surface,   5  the  boundary  layer  thickness,   u 
the  velocity  in  the  boimdary  layer,  uj  the  velocity  outside  the  boimdary 
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layer,  —  =  -pu-[_ (d.u]_ /dsj   the  pressure  gradient,  and  tq  =  -)ifSu/Syj 

the  skin  friction  at  the  surface.  Assuming  the  surface  (y  =  O)  is  made 
porous,  through  which  the  fluid  is  sucked  with  the  velocity  c,  we  have 
the  equation  of  continuity 


w  -  c  =  —  /   u  dy 
ds  J  0 


(a) 


and  the  equation  of  momentum 


punW  -  p  —  /   u'^  dy  =  TQ  +  6  — 


ds 


ds 


0^) 


where  w  is  the  velocity  of  fluid  entering  the  boundary  layer  through 
y  =  6.  The  eqixation  of  motion  reduces  to 


y=o 


(c) 


for  y  =  0. 

Now,  the  velocity  profile  in  the  boundary  layer  may  be  approximated 
by  the  Pohlhausen  polynomial 


u=Un42|^-2^+22- 


(d) 


J 


when  neither  pressure  gradient  nor  suction  exists.  For  this  profile  we 
have 


Tq  =  2H 


PU]_W  = 


1^ 


(e) 
(f) 
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To  simplify  the  calculation,  we  assume  that  the  expressions  (d) ,  (e) , 
and  (f)  still  hold  when  both  pressure  gradient  and  suction  exist.  Then 
we  have  from  (b) 


Integrating  we  have 


(g) 


d5 

630   V 

lOU  s  '^^l 

ds 

37  U16 

367  U3_  ds 

Ui"V=i-l"^l     ^^^^/   -^""'^^  (^) 


208 
where  m  =  ,  and  s,-,   is  the  initial  position  of  suction,  which  is 

367 
the  minimijm  pressure  point  in  the  present  case.  We  have  also  from  (a) 

c  = 5  (1) 

3670  ds 

If  we  substitute  (d)  and  ( i)  into  both  sides  of  the  equation  (c),  however, 
the  left  and  right  sides  become  1.003u]_/^du2_/ds  j  and  u^^dui/dsj ,  respec- 
tively.  This  contradiction  is  evidently  due  to  the  crude  assumption  of 
using  (f)  in  spite  of  the  presence  of  pressure  gradient  and  suction,  but 
we  may  overlook  the  error  because  it  is  small. 

26.  We  apply  the  calculation  to  the  symmetrical  laminar-flow  air- 
foil of  10  percent  maximum  thickness,  L.B.  2k,    set  at  zero  angle  of 
attack.   The  velocity  distribution  U]_/V  calculated  from  the  potential 

flow  of  ideal  fluids  is  used,°  the  maximum  velocity  (minimum  pressure) 
being  located  at  6k   percent  of  the  chord  from  the  leading  edge  ( sq  =  0.65t]. 

Applying  a  distributed  suction  downstream  of  the  minimum  pressure  point 
so  as  to  maintain  the  velocity  profile  in  the  boundary  layer  the  same  as 
at  this  point,  we  have  the  boundary  layer  thickness  5  and  the  required 

°It  is  ass\imed  that  the  velocity  distribution  is  not  affected  by 
the  suction.   Theoretically  n^     should  be  0  at  the  trailing  edge,  but 
the  distribution  was  somewhat  modified  so  as  to  give  uj  =.  O.85V  there. 

The  effects  of  these  assumptions  appear  to  be  too  small  to  affect  the 
result  materially. 
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suction  velocity  c  as  shown  in  figure  25.   It  is  to  "be  noted  that  s 
is  measured  along  the  surface  from  the  leading  edge,  and  R   is  the 
Reynolds  number  based  on  the  chord  length  t  and  the  velocity  of  the 
undisturbed  stream  V.   Integrating  the  area  under  the  curve  of  c,  we 

have  the  total  amoimt  of  suction  1.3vv'R  per  unit  span  of  the  two  sur- 
faces.  If  we  assume  span  =  35  m,   t  =  5  m,  V  =  200  m/s,   v  =  0.15  cm^/s, 
the  total  amount  of  suction  amounts  to   5.6  m?/a,   which  will  require  an 
exit  area  of  only  0.028  m^  when  discharged  with  the  velocity  equal  to  V, 
Integration  of  tq  gives  the  drag  coefficient  Cj)q  =  O.OOO5.  This  value 

may  be  compared  with  Cj)q  =  O.OOO3  for  the  flat  plate  with  laminar 
boundary  layer,  Cj)q  =  0.00i<-U  for  the  flat  plate  with  turbulent  boundary 

layer,  and  CpQ  =  0.0025  for  L.B.  2k   without  suction.   If  the  thickness 

of  the  airfoil  is  doubled  (20  percent  chord),  then  the  amount  of  suction 
will  be  nearly  doubled;  the  drag  is  however  almost  imchanged. 

It  should  be  noticed  again  that  the  calculation  is  based  on  the 
assumption  that  no  transition  occurs  if  the  velocity  profile  in  the 
boundary  layer  maintains  the  form  at  the  minimum  pressure.'  It  is  the 
purpose  of  the  calculation  to  show  that  extraordinarily  low  profile  drag 
may  be  expected  with  a  relatively  small  amo\ant  of  suction  under  such  a 
condition. 


PREMATURE  TRANSITION  OF  BOUNDARY  LAYER  -  EFFECT  OF  STREAM  TURBULENCE 


27.  Although  the  transition  of  the  boundary  layer  occurs  only  down- 
stream of  the  minimum  pressure,  so  far  as  flight  experiments  on  smooth 
airfoil  surfaces  are  concerned,  there  are  many  examples  of  wind  tunnel 
experiments  in  which  the  transition  moves  upstream  of  the  minimum  pres- 
sure. This  seems  to  be  due  to  the  premature  transition  ca\ised  by  the 
tiorbulence  in  the  wind  tunnel  stream.  For  example,  the  transition  on 
the  airfoil  L.B.  2k   was  found  only  downstream  of  the  minimum  pressure 
in  the  range  of  Reynolds  numbers  covered  by  the  author  's  wind  tunnel 

experiments  (the  Reynolds  number  based  on  chord  length  up  to  3  X  10°; 

see  paragraphs  12  to  1^*-) ;  as  a  result  very  low  values  of  the  profile  drag 

coefficient  C-Qq     were  observed.  The  same  airfoil,  however,  when  tested 

with  a  larger  model  of  2  m  chord  in  the  2.5  m  wind  tunnel  of  Kawasaki 
Aircraft  Company,  Gihu,  gave  the  resiilt  as  shown  in  figure  26,  in  which 

CpQ  increases  considerably  when  the  Reynolds  number  exceeds  5  X  10 

(refs.  2k   and  25).   There  is  reason  to  believe  that  the  increase  in  drag 
is  due  to  the  effect  of  stream  turbulence.   The  boundary  layer  observa- 
tion at  the  Kawasaki  wind  tunnel  shows  that  the  transition  is  found  at 

50  percent  chord  (x  =  O.5)  for  the  Reynolds  number  6  X  10°  and  moves 
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further  foirward  as  the  Reynolds  number  increases.   The  boundary  layer 
velocity  profile  observed  at  transition  has  a  form  factor  —  =  2.6 

(5*  and  9     are  the  displacement  and  momentum  thicknesses  of  the 
boundary  layer) ,  which  is  very  near  to  the  value  for  the  case  of  zero 
press-ure  gradient.   This  result  seems  to  suggest  that  the  transition  is 
at  least  not  correlated  with  the  laminar  separation  (ref .  l) . 

28.  In  order  to  verify  this  conjecture,  it  is  desirable  to  show 
that  the  transition  in  the  boimdary  layer  along  a  flat  plate  occurs 
under  the  same  condition,  because  that  transition  may  be  considered  to 
be  independent  of  the  laminar  separation.  Unfortunately,  however,  no 
flat  plate  was  measured  in  the  Kawasaki  wind  tunnel.   Therefore  we  pro- 
ceed in  a  somewhat  indirect  way.  We  assume  that  the  degree  of  stream 
turbulence  is  represented  by  the  conventional  critical  Reynolds  number 
of  the  sphere,  Rq ,  and  the  condition  of  transition  due  to  turbulence 


represented  by  the  local  Reynolds  number,  Rq 


n-^e 


at  transition  on  a 


flat  plate,  where  uj  is  the  velocity  outside  the  boundary  layer  and 
6     is  the  momentum  thickness  of  the  boundary  layer.  It  is  generally 
accepted  that  the  turbulence  in  the  wind  tunnel  stream  will  give  a 
fluctuation  of  pressure  gradient,  as  a  result  of  which  an  instantaneous 
and  intermittent  separation  will  occior.  Such  an  instantaneous  and  inter- 
mittent separation,  however,  does  not  necessarily  lead  to  the  transition 
into  t\irbulencej  for  the  transition  really  to  occur,  it  seems  probably 
necessary  that  the  Reynolds  number  Eq     which  represents  the  ratio  of 

inertia  puj^   to  viscous  stress   pvu^^/S   exceed  a  certain  critical  value. 

It  is  also  expected  that  the  critical  value  depends  on  the  degree  of  tur- 
bulence; it  must  increase  as  R^   increases.  This  is  really  shown  by 
the  experimental  data  hitherto  published,  which  are  given  in  the  following 
table  and  also  by  white  circles  (o)  in  figure  27.  The  available  data  are 
scanty,  especially  because  the  experiment  on  a  flat  plate  is  very  diffi- 
cult.  It  was  necessary  for  the  author  to  perform  a  new  experiment  (ref.  30) 
with  a  view  to  adding  one  point  in  the  range  of  high  Re. 


Rc 

Re 

Wind  Tunnel 

Reference 

i.Uo  X  io5 

0.21  X  io3 

National  Bureau  of  Standards 

26 

2.75  X  io5 

.70  X  io3 

National  Bureau  of  Standards 

26 

2.20  X  io5 

.U2  X  io3 

N.P.L.  Compressed  Air  Tunnel 

27,  28 

3.66  X  io5 

1.05  X  io3 

Aero.  Res.  Inst.  I.5  m  Tunnel 

29,  30 
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29.  Now,  we  calculate  the  value  of  Rg  at  transition  of  the  air- 
foil L.B.  2k  tested  in  the  Kawasaki  wind  timnel,  and  correlate  it  with 
the  critical  Reynolds  number  Re   of  that  tunnel.  We  analyze  similarly 

the  other  available  data,  and  summarize  the  result  in  the  following 
table.   The  values  of  the  form  factor  5*/9 ,  not  shown,  were  all  found 
in  the  range  from  2.1  to  2.7.  We  then  plot  the  data  also  in  figxire  27 
by  black  circles  (•) .  The  black  circles  are  seen  to  define  a  single 
curve  together  with  the  white  circles  already  mentioned.   This  result 
seems  to  support  the  supposition  that  the  transition  under  consideration 
is  mainly  caused  by  the  stream  turbulence,  but  not  correlated  with  the 
laminar  separation. 


Rc 

% 

Model 

Wind  Tunnel 

Reference 

1.85  X  io5 

O.i+l  X  io3 

Symmetrical    airfoil 

N.P.L.    7  ft 

31,   32 

2.10  X  io5 

.56  X  io3 

Airship  model 

M.I.T.    tI  ft 

33,  3^^ 

3.50  X  io5 

.95  X  lo3 

Airfoil  N-22 

NACA  Full-Scale 

12,  35 

3.65  X  lo5 

1.08  X  lo3 

Airfoil  L.B.   2k 

Kawasaki  2.5m 

2k,  25 

In  reference  12  (the  third  line  in  the  preceding  table) ,  the  same 
airfoil  was  examined  both  by  the  full-scale  wind  tunnel  and  by  the  flight 
tests.  We  calculate  the  form  factor  b*/9      from  these  tests  and  plot  the 
values  against  s/t  in  figure  28,  where  t  is  the  chord  length  and  s 
is  the  length  measured  along  the  surface  from  the  leading  edge.   The 
value  of  6*/0   at  transition  is  2.6  in  the  wind  tunnel,  while  it  reaches 
as  high  as  3«1  and  drops  sharply  in  the  flight  test.   The  minimum  pres- 
sure is  located  at  x-  =  O.I8,  and  the  laminar  separation' calculated  from 

the  measured  distribution  of  pressure  at  t-  =  O.36.  This  example  is  very 

interesting  because  the  cause  of  transition  is  quite  different  in  the  two 
cases  (namely,  it  is  due  to  the  stream  turbulence  in  the  wind  tunnel, 
while  it  is  related  to  laminar  separation  in  the  flight  test) ,  although 
the  positions  of  transition  are  almost  the  same. 


PREMATURE  TRANSITION  OF  BOUNDARY  LAYER  -  EFFECT  OF  SURFACE  ROUGHNESS 


30.  Up  to  this  point,  we  have  only  considered  the  case  when  the  sur- 
face  of  the  airfoil  is  smooth.   If  the  surface  is  rough,  however,  there 
is  a  possibility  that  the  transition  may  also  be  caused  prematurely  by 
surface  roughness.   So,  it  is  important  In  practice  to  estimate  the 
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approximate  order  of  magnitude  of  the  permissible  roughness  in  the  lami- 
nar boundary  layer.  Nothing  has  been  known  concerning  this  problem, 
except  a  mere  conjecture  or  fragmentary  data.   Schiller  (ref .  36)  sug- 
gested that  a  local  separation  occurs  and  hence  leads  to  transition 
when  the  Reynolds  n\imber  kuj^V  exceeds  a  certain  critical  value  ^crit' 
where  k  is  the  height  of  projection  and  Uj^  is  the  velocity  at  the 
top  of  projection.  The  exact  value  of  Rcj-it  is  ^o't  known,  but  it  is 
not  likely  to  differ  much  from  the  critical  value  of  the  Reynolds  number, 
above  which  vortices  are  shed  from  the  obstacle  of  the  same  shape  as  that 
of  the  projection  placed  in  a  uniform  stream.   The  experimental  result 
due  to  Wieselsberger  (ref.  37)  shows  that  such  a  critical  Reynolds  niomber 
is  roughly  50  for  a  circular  cylinder.  Assuming  that  the  height  of  pro- 
jection k  is  small,  and  that  the  presence  of  the  projection  in  no  way 
alters  the  character  of  the  flow,  we  have  the  shearing  stress  at  the 


surface  tq  =  pv(uij^/k).  Using  the  so-called  friction  velocity  v^  =  ^tq/p 


instead  of  Uj^,  we  have  then  kv„/v  =  ykuj^/v.  The  permissible  roughness 

is  therefore  given  by  kv^/v  =  JR^^^^,  or,  with  I^crit  ~  "^ >      ^'^■x-/'^  ~  7- 

On  the  other  hand,  according  to  Nikuradse's  experiments  on  roughened 

kv 
pipes  (ref,  38)  >  "the  critical  Reynolds  number  is  — —  =  \,   above  which 

the  roughness  projections  disturb  the  laminar  sublayer  of  the  turbulent 
boundary  layer  and  hence  increase  the  pressure  drop.   It  appears  there- 
fore that  the  permissible  roughness  is  smaller  in  the  turbulent  boundary 
layer  as  compared  with  the  laminar  boundary  layer.   This  is  confined  by 
a  British  flight  experiment  (ref.  39)  on   the  airfoil  section  of  10  feet 

chord  (Reynolds  number  1.8  X  10 7) ,  because  the  effect  of  camouflage  paint 
of  0.001  inch  thickness  increased  the  drag  by  about  6  percent  without 
moving  the  tiransition  forward.  At  any  rate,  however,  such  a  estimate  is 
nothing  but  mere  conjecture.  With  a  view  to  making  the  estimate  more 
definite,  we  performed  wind  tunnel  experiments,  although  of  small  scale, 
of  quantitative  character  (paragraphs  3I  to  32). 

31.  A  polished  aluminum  plate,  80  cm  long,  60  cm  wide,  and  3  n™  thick, 
was  held  horizontally  in  the  1.5  m  wind  tunnel  of  the  Aeronautical  Research 
Institute.   So  that  the  flow  at  entry  would  not  be  distiirbed,  the  leading 
edge  of  the  plate  was  rounded,  and  the  plate  slightly  tilted  so  that  the 
forward  stagnation  point  was  on  the  same  surface  as  that  where  the  observa- 
tion was  made.  The  tilting,  however,  was  so  slight  that  the  static  pres- 
sure was  observed  to  be  practically  uniform  along  the  plate.   The  plate 
was  roughened  by  a  wire,  which  was  stretched  across  the  flow,  in  contact 
with  the  plate.   The  diameters  k  of  the  wire  were  0.25,  0.4,  and  0.7  mm, 
respectively,  and  the  distances  x  of  the  wire  from  the  leading  edge  were 
15,  30,  ^5,  and  60  cm,  respectively.   When  the  wind  speed  V  was  low,  the 
boimdary  layer  was  laminar  all  along  the  plate,  but  from  a  certain  speed 
upward,  the  transition  to  turbulent  flow  was  observed  at  that  point  where 
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the  wire  was  placed.   Transition  was  detected  by  a  sudden  change  in  the 
value  of  the  total  pressure  G*  in  terms  of  the  dynamic  pressure  q  of 
the  undisturbed  stream,  G*  being  indicated  by  a  small  pitot  tube  with 
a  flattened  mouth  of  1  mm  external  width  and  O.3  mm  width,  which  was 
placed  in  contact  with  the  plate  at  a  point  7O  cm  behind  the  leading 
edge.  A  sample  record  of  measurements  for  k  =  0.^4-  mm  is  shown  in 
figure  29. 

When  a  flat  plate  is  placed  along  a  uniform  stream  of  velocity  V, 
the  Blasius  solution  (ref .  kO)    of  the  laminar  boundary  layer-  equation 
gives 


v^  =  0.576v(X2Ly 


■3A 


for  a  point  of  distance  x  from  the  leading  edge.  Writing  K  for  the 
critical  value  of  kv^^/v,  the  permissible  height  of  projection  k  is 

given  by 


0.576|=k(Sl)- 


-3A 


We  determine  V  from  the  wind  speed  corresponding  to  the  kink  of  the 
curve  as  shown  in  figure  29,  plot  0.5T6(k/x)   in  a  logarithmic  scale 
against  Vx/v,  and  draw  a  straight  line  of  the  slope  -3A  through  the 
points  (fig.  30) •  We  thus  obtain  K  =  I3,  which  is  far  greater  than 
the  value  K  =  7  estimated  previously. 

32,  Similar  measurements  were  also  performed  on  an  airfoil  section 
L.B.  2k.     The  model  was  of  0.8  m  span,  of  1.2  m  chord,  fitted  with  end 
plates  1.3  m  X  0.6  m,  and  set  at  zero  angle  of  attack  in  the  same  wind 
tunnel.  Wires  of  various  diameters  (k  =  0.25,  O.k,   O.7  nm)  were  attached 
parallel  to  the  span,  in  contact  with  the  surface,  at  10  percent  of  the 
chord  from  the  leading  edge  (x  =  O.l).   Transition  was  detected  by  the 
sudden  change  in  total  pressure  as  indicated  by  a  pitot  tube  with  a 
flattened  mouth  of  2.7  nm  external  width  and  O.9  mm  depth,  which  was 
placed  in  contact  with  the  surface  at  50  percent  of  the  chord  from  the 
leading  edge.  Results  of  measurements  are  shown  in  figure  31^  where 

R  =  •^,  and  t  is  the  chord  length. 

The  friction  velocity  may  be  generally  expressed  in  the  form 

v^  =  AVE"-'-' 
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where  A  is  a  function  of  s/t   (s   is  the  length  measured  along  the 
surface  from  the  forward  stagnation  point) .  We  can  calculate  A  by 
applying  either  the  Pohlhausen  approximate  solution  (ref.  4l)  or  the 
simplified  method  due  to  the  author  (refs.  15  and  16) .  The  permissible 
height  of  projection  is  then  determined  by 

A  I  =  KR-3A 

Applying  the  Pohlhausen  method  to  the  theoretically  calculated 
distribution  of  pressure,  we  get  the  values  of  function  A,  as  shown 
in  figure  32.  Since  A  is  1.23  at  the  position  of  wire  (x  =  O.l), 
the  value  1.23(k/t)   is  plotted  in  a  logarithmic  scale  against  R  in 
figure  33^  ^     being  the  Reynolds  number  corresponding  to  the  kink  of 
the  curve  as  given  in  figure  31>   I't  will  be  seen  that  although  the 
measxared  points  are  on  a  straight  line  of  the  slope  -3A^  they  give 
K  =  15,  which  is  somewhat  higher  than  the  value  found  for  the  flat  plate. 

33*  Now  we  apply  the  preceding  result  to  the  fragmentary  data 
hitherto  known  in  order  to  check  the  adequacy  of  the  estimate.  First, 
we  examine  the  results  of  wind  tunnel  experiment  on  a  symmetrical 
laminar-flow  airfoil  L.B.  27,"^  on  which  various  projections  are  attached 
at  3  percent  of  the  chord  from  the  leading  edge  (x  =  O.O3) •  The  model 
was  of  0.8  m  span,  of  1.2  m  chord,  and  set  at  zero  angle  of  attack. 
The  profile  drag  was  measured  by  the  method  similar  to  that  for  L.B.  2k 
(see  paragraphs  12  to  I3) .  The  resxilts  are  shown  in  figure  3^^  from 
which  we  find  that  the  rubber  tape  of  O.O7  mm  thickness  gives  no  effect 
over  the  range  of  Reynolds  numbers  R  covered  by  the  experiment,  while 
the  piano  wire  of  0,5  mm  diameter  gives  a  completely  turbulent  friction. 
The  effect  of  the  wire  of  0.25  inm  diameter  begins  to  appear  at 

R  =  1.3  X  10°.  Inserting  the  values  k  =  0.25  inm,  t  =  1.2  m  and 
A  =  1.95  in  the  formula  A  |-  =  KR"3A,  we  get  K  =  I5.5,  which  is  in 

good  agreement  with  the  result  in  paragraph  32.  The  value  A  =  1.95 
was  read  from  figure  32,  since  the  leading  edge  portion  of  L.B.  27 
almost  coincides  with  that  of  L.B.  2k.      Profile  drag  coefficients  of 
both  airfoils  are  also  the  same  over  the  range  of  Reynolds  numbers 
examined . 

3^.  In  connection  with  the  determination  of  boundary  layer  transi- 
tion on  airfoils  in  the  NACA  full-scale  wind  tunnel  (ref.  k2) ,   an  aux- 
iliary measiirement  has  been  reported,  in  which  the  effect  was  examined 

7 

'L.B.  27  has  a  maximum  thickness  of  10  percent  of  the  chord  at 

6D  percent  of  the  chord  from  the  leading  edge.  See  figure  5  of  paragraph  6. 
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of  rubber  tapes  attached  at  5  percent  of  the  chord  (x  =  O.O5)  from  the 
leading  edge  of  the  model.  The  airfoil  section  was  NACA  0012,  the  chord 

was  72  inches,  and  the  Reynolds  number  R  was  U.I8  X  10°.  No  effect 
was  foimd  when  the  tape  was  O.OO3  inch  thick,  some  effect  began  to 
appear  when  it  was  O.OO6  inch  thick,  and  the  transition  moved  right  to 
that  position  where  the  tape  of  O.OO9  inch  thickness  was  attached. 
Assuming  A  =  1.6  and  K  =  15,  we  estimate  from  the  preceding  formula 
the  value  O.OO7  inch  for  the  permissible  thickness,  which  seems  to  agree 
well  with  the  observation. 

If  we  further  ass\me  that  the  value  A  =  1.6  is  also  applicable 
to  the  case  of  the  British  flight  experiment  mentioned  in  paragraph  30^ 
we  find  0.004  inch  for  the  critical  height  for  transition  with 

t  =  10  feet,  R  =  1.8  X  10'^,  and  K  =  I5.   On  the  other  hand,  we  esti- 

kv„   , 
mate  the  permissible  limit  in  turbulent  boundary  layer  by  — —  =  ^^ 


which  may  be  written  in  the  form 


k 

t 


^c^  u-^  R 


PI     P 
by  the  relation  v^  =  ±-  c^uj^ ,   where  uj^     is  the  local  wind  speed  and 

c-f  is  the  coefficient  of  local  skin  friction.   If  we  assume  u-]_  =  1.2V 

and  Cf   =  0.003  (which  is  equal  to  the  coefficient  of  mean  skin  friction 

for  a  flat  plate  at  R  =  1.8  X  10 "7),  we  obtain  O.OOO6  inch  for  the  per- 
missible roughness  thickness.  Since  the  thickness  of  the,  camoiiflage 
paint  is  reported  to  be  about  0.001  inch,  it  may  be  concluded  that  the 
paint  increases  only  the  friction  in  turbulent  boundary  layer,  without 
affecting,  however,  the  transition  to  turbulent  flow.  This  is  in  good 
agreement  with  the  experimental  results. 
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APPENDIX 


TRANSITION  CAUSED  BY  LAMINAR  SEPARATION 

35.  As  is  well-known,  the  phenomenon  of  sudden  decrease  in  drag  of 

a  sphere  at  a  certain  value  of  the  Reynolds  number  R  =  —  (V  is  the 

speed  of  undist\irbed  stream  and  d  is  the  diameter  of  sphere)  is 
explained  by  supposing  that  the  boundary  layer  separates  while  it  is 
laminar  when  R  is  low,  but  it  separates  after  transition  to  turbulence 
when  R  is  high,  thus  resulting  in  diminishing  the  so-called  dead  water 
region.   Probably  the  transformation  from  laminar  separation  to  turbulent 
separation  may  proceed  as  follows: 

When  the  laminar  boundary  layer  separates  from  the  surface,  the 
detached  layer  remains  also  laminar  at  first,  but  it  is  so  unstable  that 
it  becomes  turbulent  at  a  short  distance.  This  transition  from  laminar 
to  turbulent  flow  is  considered  to  occur  when  the  local  Reynolds  number 
based  on  the  width  of  the  detached  layer  and  the  velocity  outside  the 
layer  exceeds  a  certain  value,  so  that  the  transition  moves  upstream 
toward  the  separation  point  as  R  increases.  When  the  transition 
approaches  sufficiently  near  the  separation  point,  it  becomes  possible 
for  the  detached  layer  to  come  back  again  to  the  downstream  surface, 
because  the  tiirbulence  produced  will  drive  the  flow  forward.   The  layer 
reattaches  to  the  siirface  as  a  tiirbulent  layer,  and  accordingly  the  drag 
coefficient  begins  to  decrease.  The  distance  between  the  separation  and 
the  first  turbulent  boundary  layer  decreases  as  R  increases,  and  finally 
the  fully  developed  turbulent  boundary  layer  commences  just  downstream 
of  the  separation  point.   The  drag  coefficient  then  ceases  to  decrease. 

36.  Now,  in  order  that  the  separated  layer  reattach  to  the  surface, 

U-10 

it  seems  necessary  for  the  local  Reynolds  n\imber  Rg  =  — =—  at  separation 

to  exceed  a  certain  critical  value,  where  u-^  is  the  velocity  outside 

the  boundary  layer,  and  0   is  the  momentum  thickness  of  the  boundary 
layer.  This  may  be  explained  as  follows:  According  to  the  laminar 
boundary  layer  theory,  the  separation  occiirs  when  the  quantity 

_02dui^^'^^    ui 
V  ds    ds    •     0 

exceeds  a  certain  value,  suggesting  that  the  pressure  rise  (dp/ds)0 
becomes  too  large  in  proportion  to  the  shearing  stress  at  the  sur- 
face pv(ui/0).  Assimiing  analogically  that  the  separated  layer  leaves 
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35 


the  siorface  when  the  pressure  rise  becomes  too  large  in  proportion  to 
the  momentum  pu-i  ,  we  then  find  that  Rq   at  the  separation  point  must 

exceed  a  certain  critical  value  in  order  that  the  separated  layer 
reattach  to  the  surface. 

37»  In  order  to  determine  the  critical  value  of  Rq ,  we  consider 

in  detail  the  condition  where  the  coefficient  of  sphere  drag  begins  to 
decrease.   This  condition  corresponds  to  the  point  B  of  the  curve  of 
figure  35^  which  represents  an  idealized  variation  of  the  drag  coeffi- 


;ient  C^p  drag  -J-  (pv2/2)  fnd^/l^j 
?/qU  di 


or  the  pressure  difference  coeff i- 


with 


cient  Ap/q  =  difference  of  pressures  at  the  forward  stagnation  point 

and  the  point  corresponding  to  the  central  angle  157.5°  -7  (pV  /2) 

the  Reynolds  number  R.  Within  the  range  AB,  the  pressure  distribution 
aroimd  the  sphere  is  approximately  independent  of  R;  the  typical  example 
may  be  found  from  the  experiments  due  to  Fage  (ref .  hs) .      Fortunately, 
the  boundary  layer  calculation  has  also  been  performed  for  that  distri- 
bution of  pressure  by  Tomotika  and  Imai  (ref.  hk) ,   so  that  the  local 
Reynolds  number  Rg   is  given  by 

Re  =  o.i+o\Ir 

at  the  separation  point.  Althoiigh  the  calculation  has  originally  been 

made  for  a  particular  Reynolds  number,  R  =  1.57  X  10^^  the  preceding 
relation  may  be  applied  for  any  value  of  R  in  the  range  of  AB.   Putting 
the  value  of  R  at  B,  and  writing 

^crit  =  °-^°^ 

we  have  Rg      as  the  critical  value  of  Rg   above  which  the  separated 

laminar  layer  reattaches  to  the  surface.  Conventionally  the  Reynolds 
number  Rq   corresponding  to  Cj,  =  0.3   (or  Ap/q  =  1.22)   has  been  used, 

instead  of  R^,  for  representing  the  degree  of  stream  turbulence,  but  it 
is  not  so  difficult  to  estimate  the  value  of  Rg  from  the  measured  c\irve 
of  Cj)   (or  Ap/q)   against  R.  For  example,  we  have  from  the  experiments 
of  towing  spheres  in  the  free  atmosphere  (ref.  35) 

Rg  =  3.6  X  10^     Rq  =  3.85  X  105     Rg   _   =  2l|0 


36 
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Since  these  values  refer  to  the  case  of  very  low  turbulence,  it  will  be 

seen  that  Rq     =  2^*0  represents  the  highest  possible  value.   It  is 
crit 

also  to  be  noted  that  -=■  =0,9^  in  this  example  and  that  almost  the 

same  value  has  been  obtained  by  the  author's  experiments  on  spheres  of 
various  diameters  (ref .  29) . 

38.  If  the  stream  turbiilence  is  not  low  so  that  Bq      is  less  than 
3.85  X  10^,  then  Eq  will  be  less  than  2^0.  Assuming  the  ratio  Rg/Rc 

to  be  constant,  we  can  estimate  the  corresponding  value  by 


Rfl     =  2l4-0| 
^crit 


R^ 


lis. 85  X  10^ 


On  the  other  hand,  we  can  also  estimate  the  value  of  %ri.it  directly 

from  the  boundary  layer  measurements.   The  results  of  the  analysis  for  a 
sphere  as  well  as  circular  and  elliptic  cylinders  are  summarized  in  the 
following  table,  where  Rq      is  the  critical  value  Rq  estimated 

from  Rp  by  the  preceding  formula,  Rq     is  the  value  of  Ra   observed 

sep 

at  the  separation  point  when  the  boundary  layer  really  separates  while 

it  is  laminar,  and  Ba  is  the  value  of  Ra   observed  at  the  calcu- 

" trans 

lated  laminar  separation  point  when  the  boundary  layer  separates  after 

transition.  The  fact  that  Ea      .,   lies  between  R^     and  Rfl 

^crit  ^sep       ^trans 

seems  to  suggest  the  adequacy  of  the  preceding  consideration. 


Body 

^ 

^crit 

^sep 

■^  trans 

Reference 

Sphere 

Circular  cylinder 

Elliptical  cylinder 

2.5  X  io5 
1.5  X  lo5 
2.7  X  io5 

190 
150 
200 

160 
lUO 
160 

220 

225 
1+00 

^3 

1^5 

^6,  1+7 

39 •  We  now  proceed  to  apply  our  result  to  inteiTpreting  the  effect  of 

Reynolds  nijmber  on  maximum  lift  of  airfoils.   For  the  angle  of  attack 
near  the  stall,  the  flow  separates  shortly  downstream  of  the  leading  edge 
while  the  boundary  layer  is  laminar.   If  the  flow  fails  to  reattach  to 
the  surface  as  a  turbulent  layer,  the  maximum  lift  coefficient  Ct^   of 
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the  airfoil  will  be  almost  independent  of  the  Reynolds  nimiher  R.  At 
the  separation  point,  similar  to  the  case  of  a  sphere,  the  relation  of 
the  form 


Rfl  =  klfR 


holds,  where  R  =  —  is  the  Reynolds  number  referred  to  chord  length  t, 

and  k  is  a  constant  depending  on  the  shape  of  airfoil  and  the  value  of 

Ct_   .   If  R  is  low  so  that  Rq   is  less  than  the  critical  value  R^ 
^-max  ^  ^crit 

then  Ct     will  be  independent  of  R.  Assuming  the  same  value  of 

Rfl      for  the  sphere  as  for  the  airfoil,  we  obtain 
^crit  ^  ' 

for  the  Reynolds  number  above  which  Ct     begins  to  increase  with  the 

Reynolds  number.   Therefore,  the  ratio  of  the  Reynolds  number  corresponding 
to  a  certain  value  of  Ct     of  an  airfoil  and  the  critical  Reynolds 

niomber  of  sphere  in  the  same  stream,  Rg  or  Rq ,  becomes  independent  of 
the  stream  turbxilence.  Denoting  the  values  for  a  reference  tunnel  with 
asterisk,  we  have 

5I  =  ^ 
R   Re 

which  in  tiorn  means  that  the  ratio  of  Reynolds  numbers  corresponding  to 
a  certain  value  of  C-r     is  equal  to  the  ratio  of  critical  Reynolds 

numbers  of  a  sphere.   This  is  useful  for  comparing  the  values  of  ^Lnjoy 

obtained  in  two  different  wind  tunnels.  Considering  the  reference  condi- 
tion to  be  the  free  flight  in  the  atmosphere,  we  find 

R*  =  R  X  -2_ 
RC 


°It  is  to  be  noted  that  R  is  referred  to  t,  while  Rg  and  Bq 
are  referred  to  d. 
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as  the  free  flight  Reynolds  number  which  will  give  the  same  value  of 
Cl^    as  that  observed  in  a  wind  tunnel.  This  is  just  what  is  called 

the  effective  Reynolds  number.   Strictly  speaking,  such  an  argument  as 
mentioned  before  should  apply  only  to  the  Reynolds  numbers  near  the 
critical  value,  but  there  are  many  experimental  evidences  showing  the 
usefulness  of  the  concept  of  effective  Reynolds  number  for  most  practi- 
cal purposes,  as  far  as  the  commonly  used  airfoils  and  range  of  Reynolds 
numbers  of  both  wind  tunnel  and  free  flight  are  concerned. 

kO.   Finally,  we  consider  a  more  quantitative  example  to  show  the 

adequacy  of  the  preceding  argument.   In  figure  36,  Cj     for  various 

■Tnax 

NACA  symmetrical  airfoils  are  plotted  against  the  effective  Reynolds 

number  R*,  the  experimental  data  being  taken  from  the  results  of  the 

NACA  variable-density  wind  tvmnel  (ref.  U8) .  As  already  mentioned,  up 

to  a  certain  value  of  R*,  Ct     is  almost  independent  of  R*.  This 

corresponds  to  the  condition  in  which  the  laminar  separation  just  behind 
the  leading  edge  fails  to  reattach  to  the  surface,  resulting  in  a  con- 
siderable dead  water  region  above  the  airfoil  surface.  The  value  of 
Cr     is  approximately  O.9,  irrespective  of  the  thickness;  it  is  almost 

equal  to  the  value  for  a  flat  plate  of  vanishing  thickness.  Theory  of 
discontinous  flow,  when  applied  to  the  flat  plate,  seems  to  give  a  lift 
coefficient  close  to  O.9  (ref.  ^4^9)  •  We  idealize,  therefore,  the  experi- 
mental curve  as  shown  by  dotted  lines  in  figure  36.  Then,  the  point 
where  the  dotted  line  meets  the  line  Ct    =0.9  will  be  considered  to 

correspond  to  Rp     =  2^0.   In  order  to  determine  this  point,  we  calcu- 
^crit 

late  the  value  of  R0/\/R  at  the  laminar  separation  point  for  a  lift 

coefficient  Cl  =  O.9,  and  the  value  of  R  which  gives  Rq   =  2l+0.  We 

first  calculate  the  pressure  distribution  by  the  formula  of  paragraph  3 

for  the  two-dimensional  potential  flow  around  the  airfoil  section. 

Although  the  formula  may  be  applied  to  any  arbitrary  airfoil  section, 

we  have  determined  the  pressure  distribution  only  for  the  airfoil 


T  =  0.287e^x(l  -  x)(5  -  ^x) 


in  order  to  simplify  the  calculation,  because  no  great  exactitude  is 
required  in  the  present  problem,   x  is  the  coordinate  along  the  chord, 
X  =  0  and  x  =  1  corresponding  to  the  leading  and  trailing  edges, 
respectively,  T   is  the  half -thickness,  and  e  is  the  maximum  thick- 
ness in  terms  of  chord  length.   The  airfoil  represented  by  the  preceding 
expression  coincides  with  siifficient  accuracy  with  the  true  NACA  symmetri- 
cal airfoil  except  near  the  trailing  edge.  The  value  of  R^/^  at  the 

laminar  separation  point  was  then  determined  for  the  calculated  pressure 
distribution  by  applying  the  approximate  method  due  to  the  author  (refs.  15 
and  16) . 
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The  critical  values  R    .,   th\is  calculated  are  shown  hy  a  solid 

crit  '' 

line  in  figure  37,  while  the  corresponding  values  taken  from  figure  36 

are  shown  by  white  circles.   The  agreement  is  fairly  good,  and  especially 

satisfactory  when  the  thickness  of  the  airfoil  is  small.   In  general,  the 

thickness  of  the  boundary  layer  near  the  trailing  edge  increases  as  C-^ 

increases.   If  Cj^  is  further  increased,  however,  a  laminar  separation 

suddenly  occiirs  near  the  leading  edge  when  the  thickness  is  small,  while 
the  trailing-edge  turbulent  separation  moves  a  considerable  extent  for- 
ward before  the  leading-edge  laminar  separation  occurs  when  the  thickness 
is  moderate.  Therefore,  the  asstnnption  of  the  analysis  is  more  satis- 
factorily realized  in  the  case  of  small  thickness,  thus  bringing  the 
calculated  and  observed  values  in  close  agreement. 

In  conclusion,  the  author  wishes  to  acknowledge  his  indebtedness 
for  the  assistance  given  by  Messrs.  C.  Noda,  S.  Mituisi,  I.  Shinra, 
S.  Asaka,  R.  Hama,  and  K.  Takeda. 
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TABLE  1,-  ORDINATES  FOR  SYMMETRICAL  SECTION.   VALUES  OF  T/e 


X 

I 

J 

K 

L 

M 

N 

0 

0 

0 

0 

0 

0 

0 

,003 

.oi+65 

.0563 

.0571+ 

.O58I+ 

.0606 

.0627 

.006 

.0661 

.0793 

.0808 

.0823 

.0851+ 

.0886 

.0125 

.0962 

.1135 

.1158 

.1180 

.1226 

.1275 

.025 

.137^^ 

.1-589 

.1620 

.1652 

.1720 

.1791+ 

.05' 

.1963 

.2208 

.2252 

.2298 

.2396 

.2508 

.075 

.24llj- 

.2663 

.2716 

.2771 

.2892 

.3031 

.10 

.2789 

.3029 

.3089 

.3151 

.3288 

.31+1+7 

.15 

.339^ 

.3599 

.3667 

.3738 

.3893 

.1+072 

.20 

.3865 

.1+026 

.1+097 

.1+170 

.1+328 

.1+^3 

.25 

.1+236 

.1+353 

.1+1+22 

.1+1+92 

.1+638 

.1+788 

.30 

.i+523 

.1+601 

.1+663 

.1+721+ 

.1+81+1+ 

.1+91+9 

•  35 

.^737 

.1+783 

.1+833 

.1+881 

.1+962 

.5000 

.i^O 

.U885 

.1+906 

.1+9^2 

.1+971 

.5000 

.1+953 

.45 

.i+972 

.1+977 

.1+991+ 

.5000 

.1+91+8 

.1+818 

.50 

.5000 

.5000 

.1+99U 

.1+931 

.1+797 

.k6oh 

.55 

.4970 

.1+91+1+ 

.1+873 

.1+781 

.1+558 

.i+320 

.60 

.1+871 

.1+778 

.1^656 

.1+522 

.1+21+1 

.3975 

.65 

.i+691 

.1+509 

.1+31+2 

.1+175 

.3856 

.3578 

.70 

.1+1+18 

.1+139 

.3939 

.3751 

.31+13 

.3138 

.75 

.1+038 

.3675 

.31+55 

.3258 

.2921 

.2661+ 

.80 

.3538 

.3121 

.2899 

.2706 

.2393 

.2165 

.85 

.2908 

.2i+8l 

.2277 

.2105 

.1836 

.1650 

.90 

.2133 

.1762 

.1598 

.11+64 

.1262 

.1128 

.95 

.1201 

.0966 

.0870 

.0793 

.0680 

.0609 

1.00 

.0100 

.0100 

.0100 

.0100 

.0100 

.0100 
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TABLE  2.-  AUXILIARY  FUNCTION  Bg  ASSOCIATED  WITH  THE 
PRESSURE  DISTRIBUTION  OF  SYMMETRICAL  SECTIONS 


X 

I 

J 

K 

L 

M 

N 

0.0125 

0.59 

1.10 

1.13 

1.16 

1.20 

1.21+ 

.025 

.66 

1.06 

1.09 

1.13 

1.20 

1.26 

.05 

.75 

l.Oil 

1.07 

1.11 

1.20 

■  1.29 

.10 

.86 

1.03 

1.07 

1.11 

1.21 

1.33 

.20 

.96 

1.03 

1.08 

1.12 

1.22 

1.36 

.30 

1.02 

1.0^+ 

1.09 

I.1I+ 

1.23 

1.33 

.ko 

1.05 

1.05 

1.10 

1.15 

1.23 

1.22 

.50 

1.08 

1.17 

1.19 

1.19 

l.lU 

l.OU 

.60 

1.12 

1.18 

1.11 

1.05 

.92 

.80 

.70 

1.09 

.98 

.88 

.79 

.63 

.51 

.80 

.89 

.65 

.52 

.k2 

.28 

.19 

.90 

.33 

.06 

-.03 

-.10 

-.18 

-.22 

.95 

-.32 

-M 

-.50 

-.52 

-.53 

-.52 

.975 

-.95 

-.93 

-.91 

-.89 

-.81+ 

-.78 
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TABLE  3.-  AUXILIARY  FUNCTION  Be  ASSOCIATED  WITH  THE 
PRESSURE  DISTRIBUTION  OF  SYMMETRICAL  SECTIONS 


X 

I 

J 

K 

L 

M 

N 

0.0125 

-k.ok 

-1^.1+0 

-i^.88 

-U.90 

-5.20 

-5.30 

.025 

-2.86 

-3.08 

-3.20 

-3.22 

-3.37 

-3.55 

.05 

-2.0U 

-2.08 

-2.13 

-2.17 

-2.28 

-2.39 

.10 

-1.39 

-1.33 

-1.35 

-1.37 

-I.U2 

-l.i;8 

.20 

-.83 

-.7h 

-.7^+ 

-.7h 

-.Ih 

-.70 

.30 

-.50 

-M 

-.1+1 

-.39 

-.31 

-.21 

.Uo 

-.2k 

-.19 

-.17 

-.12 

.00 

.18 

.50 

.00 

.00 

.11 

.22 

.39 

.^ 

.60 

.27 

.kk 

.53 

.60 

.70 

.75 

.70 

.65 

.83 

.89 

.92 

.9k 

.91 

.80 

1.11 

1.19 

1.18 

1.15 

1.09 

1.02 

.90 

1.70 

1.52 

l.i+l 

1.31 

1.16 

l.Oi^ 

.95 

2.03 

1.66 

1.50 

1.36 

1,16 

1.03 

.975 

2.20 

1.73 

1.^ 

1.38 

1.15 

1.02 
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TABLE   h.-  ORDINATES  FOR  MEAN  CAMBER  LINES.      VALUES   OF     M/f 


X 

Do 

Dl 

D3 

1^5 

Deo 

E1+ 

^3 

0 

0 

0 

0 

0 

0 

0 

0 

.003 

.0295 

.0282 

.0292 

.0300 

.0371 

.0313 

.0327 

.006 

.0529 

.0509 

.0526 

.051+0 

.0666 

.O56I+ 

.0589 

.0125 

.0969 

.0911+ 

.0969 

.0995 

.1220 

.1039 

.1081 

.025 

.1687 

.1636 

.1695 

.17^0 

.2119 

.1811+ 

.1887 

.050 

.286U 

.2797 

.2897 

.2973 

.3585 

.3093 

.3212 

.075 

.3843 

.3772 

.3907 

.1+006 

.1+787 

.1+161 

Mi^ 

.10 

.1^690 

.1+620 

.1+781+ 

.1+902 

.5808 

.5081+ 

.5266 

.15 

.6098 

.601+5 

.6252 

.6396 

.lkk-3 

.6615 

.6831 

.20 

.7219 

.7191 

.7i+2l+ 

.7579 

.8638 

.7812 

.801+1 

.25 

.8113 

.8112 

.8353 

.8507 

.9451 

.8733 

.8952 

.30 

.8813 

.8838 

.9067 

.9207 

.9898 

.9387 

.9586 

.35 

.93^+1 

.9381+ 

.9579 

.9688 

.9963 

.9826 

.9939 

.i+0 

.9710 

.9758 

.9892 

.9953 

.9589 

.9998 

.9973 

.h^ 

.9928 

.9962 

1.0000 

.9988 

.8625 

.9872 

.9679 

.50 

1.0000 

.9988 

.9881 

.9761 

.6300 

.9387 

.9085 

.55 

.9928 

.9813 

.9^+72 

.9178 

.3885 

.8561 

,8238 

.60 

.9710 

.9^14 

.8735 

.8192 

.261+6 

.7^+73 

.7201+ 

.65 

.93i^l 

.8794 

.770I1 

.6903 

.1811 

.6259 

.6056 

.70 

.8813 

.7878 

.61+59 

.5^7^^ 

.1207 

.1+918 

.1+871 

.75 

.8113 

.6898 

.5108 

.1+078 

.0772 

.3682 

.3721 

.80 

.7219 

.5673 

.3766 

.281+7 

.01+59 

.2590 

.2669 

.85 

.6098 

.1+305 

.2535 

.181+6 

.021+2 

.1686 

.1761 

.90 

.U690 

.281+5 

.11+93 

.1071 

.0101 

.0973 

.1020 

.95 

.2861+ 

.1366 

.0651 

.01+72 

.0021+ 

.01+25 

.01+1+2 

1.00 

0 

0 

0 

0 

0 

0     • 

0 

X     for     M/f  =   1 

.500 

1 

.1+82 

.1+50 

M3 

.333 

.1+06 

.381 

48 


NACA  TM  1351 


TABLE  5.-  AUXILIARY  FUNCTION  Ac  ASSOCIATED  WITH  THE 
PRESSURE  DISTRIBUTION  OF  MEAN  CAMBER  LINES 


X 

Do 

El 

D3 

I>5 

D» 

El, 

^3 

0.0125 

h.^3 

-0.53 

-2.50 

-3.17 

-5.i+9 

-3.51 

-3.62 

.025 

i^.53 

.88 

-M 

-.91 

-2.17 

-1.07 

-1.08 

.05 

1^.53 

1.89 

.99 

.71 

.22 

.67 

.1^ 

.10 

^.53 

2.61 

2.03 

1.88 

1.94 

1.93 

2.05 

.20 

i^.53 

3.15 

2.80 

2.75 

3.20 

2.85 

3.01 

.30 

U.53 

3.i+0 

3.17 

3.15 

3.79 

3.29 

3.i^6 

.ho 

1^.53 

3.56 

3.^0 

3.i^l 

1^.18 

3.57 

3.61 

.50 

i^.53 

3.68 

3.57 

3.61 

k.h6 

3.50 

3.10 

.60 

i+.53 

3.61 

3.21 

2.95 

-1.03 

2.35 

2.li+ 

.70 

U.53 

3.19 

2.07 

1.31 

-.83 

1.00 

I.OI+ 

.80 

i^.53 

2.J+3 

.75 

.05 

-.63 

.03 

.17 

.90 

i^.53 

1.3i^ 

-.05 

-.26 

-.k2 

-.26 

-.22 

.90 

4.53 

.68 

-.15 

-.20 

-.29 

-.21 

-.21 

.975 

i^.53 

.33 

-.08 

-.li^ 

-.20 

-.15 

-.15 
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TABLE  6.-  AUXILIARY  FUNCTION  Ag  ASSOCIATED  WITH  THE 
PRESSURE  DISTRIBUTION  OF  MEAN  CAMBER  LINES 


X 

EO 

Dl 

»3 

^5 

Deo 

E4 

^3 

0.0125 

-6.31 

-6.13 

-6.35 

-6.52 

-7.92 

-6.79 

-7.07 

.025 

-5.28 

-5.18 

-5.37 

-5.50 

-6.61 

-5.73 

-5.94 

.05 

-4.25 

-4.21 

-4.36 

-4.46 

-5.25 

-4.63 

-4.79 

.10 

-3.17 

-3.19 

-3.29 

-3.36 

-3.78 

-3.46 

-3.55 

.20 

-1.99 

-2.05 

-2.09 

-2.10 

-2.00 

-2.11 

-2  ..11 

.30 

-1.21 

-1.27 

-1.22 

-1.18 

-.52 

1.09 

-.99 

.40 

-.58 

-.58 

-.43 

-.31 

1.26 

-.07 

.26 

.50 

0 

.14 

.50 

.76 

00 

1.33 

1.46 

.60 

.58 

1.02 

1.79 

2.32 

2.00 

2.38 

2.21 

.70 

1.21 

1.90 

2.63 

2.87 

1.02 

2.30 

2.36 

.80 

1.99 

2.61 

2.60 

2.24 

.52 

2.00 

1.97 

.90 

3.17 

2.97 

1.88 

1.35 

.21 

1.25 

1.31 

.95 

4.25 

2.91 

1.63 

1.05 

.10 

.91 

1.01 

.975 

5.28 

2.75 

1.59 

.94 

.05 

.82 

.88 
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TABLE   7.-  PRESSURE  DIFFERENCE  DISTRIBUTION     G     FOR  MEAN  CAMBER  LINES 


X 

Do 

El 

^3 

1^5 

D„> 

E4 

^3 

0.0125 

i^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

4.93 

.025 

^^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

4.93 

.05 

^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

4.93 

.10 

h.^3 

4.22 

4.34 

4.47 

5.72 

4.70 

4.93 

.20 

i^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

^+.93 

.30 

i^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

^+.93 

.ko 

i^.53 

4.22 

4.34 

4.47 

5.72 

4.70 

4.78 

.50 

1^.53 

4.22 

4.34 

4.47 

5.72 

4.42 

4.06 

.60 

i^.53 

4.05 

3.84 

3.65 

0 

3.10 

2.92 

.70 

1^.53 

3.54 

2.57 

1.87 

0 

1.60 

1.67 

.80 

1^.53 

2.70 

1.14 

.48 

0 

M 

.65 

.90 

^.53 

1.52 

.20 

.03 

0 

.05 

.11 

.95 

^.53 

.80 

.03 

.00 

0 

.00 

.02 

.975 

1^.53 

.41 

.00 

.00 

0 

.00 

.00 
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